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ABSTRACT 



The Army RAH-66 Comanche Helicopter is currently undergoing developmental flight 
testing. The empennage of the aircraft is experiencing buffeting where the horizontal and 
vertical tail vibrate at resonant frequencies. These high buffet loads are manifested in 
higher than anticipated fitting loads, particularly on the tail, and vibrations in the crew 
stations and at the nose cone where the targeting sensors are located. Significant effort 
has been devoted to identifying the sources of excitation and the nature of the structural 
response. This thesis determines the location and magnitude of empennage vibratory 
airloads. Because the nature of the excitation is a random function, spectral analysis is 
used. To obtain the loads, a three step process was utilized. First, from aircraft 
differential pressure transducers and accelerometers, the spectral content of the response 
and excitation was determined Then, using a NASTRAN model modified to replicate the 
flight test aircraft, frequency response functions were determined between selected points 
on the aircraft's tail and the accelerometers. Finally, using this information, a solution was 
obtained for the vibratory airloads. Having provided information on the nature of the 
driving forces, structural modifications can be made that move the natural frequencies 
away from the frequencies of the applied airloads. 
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I. 



INTRODUCTION 



A. GENERAL 

The RAH-66 Comanche is the United States Army’s newest light/attack armed 
reconnaissance helicopter. It is designed to be a combat multiplier on the joint and 
combined arms battlefields of the 21 st century. The aircraft will be fielded to both corps 
and division cavalry and attack battalions/squadrons with the assigned missions of armed 
reconnaissance and attack, both with an embedded air to air capability. Mission profiles 
include day, night, reduced visibility, from ships, and in desert and non-desert 
environments. The aircraft is shown in Figure 1.1. 

Comanche will contribute to combat effectiveness by its ability to fight for 
information and disrupt the enemy's tactical intelligence collection effort in the counter 
reconnaissance role. Key attributes of the aircraft include a five-blade bearingless main 
rotor , FANTAIL anti-torque system, composite materials, Low Observables (LO), an 
advanced pilotage and targeting system, and an integrated digital mission equipment 
package (MEP). These digital systems will allow it to serve as a forward data fusion 
center and provide near real time information to commanders at all levels. 

Designed for two crew members, it is capable of being operated by a single crew 
member from either station. During missions, the pilot-vehicle-interface (PVT) is designed 
to support low workload, high situational awareness, human error resistance, and 
maximum mission productivity. [Ref. 1] 
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Figure 1.1: US Army RAH-66 Comanche Helicopter. 



The Boeing Defense and Space Group's Helicopter Division of Philadelphia, 
Pennsylvania and United Technologies’ Sikorsky Aircraft Division of Stratford, 
Connecticut were awarded the demonstration/validation (DEM/Val) phase of the contract 
for the Comanche. These contractors divided the aircraft into two sections for design and 
fabrication. Boeing has responsibility for the tail section of the aircraft. Sikorsky has 
responsibility for the forward portion of the aircraft fuselage, to include the integration of 
both sections. 

Currently, the first prototype is undergoing extensive flight testing at the Sikorsky 
test flight facility in West Palm Beach, Florida. Focused on "envelope expansion", the 
aircraft had flown approximately 91 hours, reached 171 knots in forward, level flight, 
achieved bank angles of 45°, and sideward flight of 70 knots as of May, 1998. Prototype 
#2 is scheduled to begin flight testing in 1999. 
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B. PROBLEM DEFINITION 



During flight testing the empennage of the aircraft experienced buffeting where the 
horizontal and vertical tail vibrated excessively in certain flight regimes. Observable from 
both the chase aircraft and on flight data records for strain gauges and accelerometers, 
these high buffet loads created higher than anticipated fitting loads, particularly in the tail 
region, and resulted in vibrations in the crew stations and at the nose cone where the 
targeting sensors are located. Significant effort has been devoted to identifying the 
sources of excitation and the nature of the structural response. One such method is shown 
in Figure 1.2. Here, strips of material were placed on the empennage to allow for better 
visualization of both the airflow and vibration of the structure. 




Figure 1.2: Comanche Prototype Test Flight. 
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Boeing and Sikorsky have also conducted several wind tunnel experiments to 
better quantify the flow coming from the aircraft fuselage. The first was completed in 
April 1997 at NASA Langley. Two follow-on tests were executed at the United 
Technologies’ Research Center in Hartford, Connecticut and then Boeing Defense and 
Space Group’s Helicopter Division in Philadelphia, Pennsylvania. These tests consisted of 
changing aerodynamic surfaces to modify the flow impinging on the tail from the fuselage. 
Additionally, the flight test aircraft had several structural modifications to stiffen the tail 
section; notably a bracket to secure the horizontal and vertical tail, use of a rigid flap on 
the right trailing edge of the horizontal stabilizer, and the adding of tip weights to the 
horizontal stabilizer. 

This research determines the location and magnitude of empennage loads caused 
by the random excitation. Because the nature of the excitation is assumed to be a random 
process, spectral analysis is used. To obtain the loads, a three step process is utilized. 
First, from aircraft differential pressure transducers and accelerometers, the spectral 
content of the excitation and response is determined. Then, using a NASTRAN model 
modified to replicate the flight test aircraft, frequency response functions are obtained 
between selected points on the aircraft's tail and the accelerometers. Finally, using this 
information, estimates of the airloads are calculated. Having provided information on the 
nature of the excitation, structural modifications can then be explored to reduce the 

response of the airframe to the excitation. Potential methods for accomplishing this may 

* 

be based on shifting airframe natural frequencies, increasing structural damping, or adding 
vibration control devices. 
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II. HELICOPTER VIBRATIONS 



A. GENERAL 

As with all helicopters, the RAH-66 experiences vibrations. These forces can 
significantly impact crew effectiveness, imaging and targeting performance and airframe 
and equipment reliability and maintainability. With the use of complex digital avionics 
systems, there is a greater emphasis on vibration reduction. As such, the analysis of 
dynamic loads and modification of structures to create a low vibration environment is 
essential. The first step in reducing vibrations is to identify the source of excitation. Once 
that is known, the aircraft can be modified so that the structural frequencies are offset 
from the driving frequencies. 

Helicopter vibrations fall into three categories. First, there are vibrations due to 
rotor excitation, which are at frequencies that are integral multiples of the rotor’s 
rotational speed. Second, vibrations due to random aerodynamic excitation where the 
frequency being observed is a natural frequency of the structure being excited. And 
finally, self-excited vibrations, such as flutter and ground resonance. [Ref 2] 

1 . Vibrations at Integral Multiples of Rotor Speed 
The primary source of helicopter fuselage vibrations is the rotor system. Here, 
harmonics of aerodynamic loads on the blade give rise to vibratory response of the blade. 
Since the blade is restrained at the root, the blade responses result in root shears, which 
feed from the rotor head into the fuselage as vibratory shears and moments. As the forces 
go from the rotating to the fixed fuselage system, the rotor system tends to act as a filter. 
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For an N blades rotor system, the frequencies which will filter through are those at N and 
2N/rev. Since the amplitudes of the lower harmonics of blade loading are greater than the 
higher, the N/rev vibration of the fuselage is most critical. It can be shown that N/rev 
fuselage vibrations in the fixed system are the result of the N-l, N, and N+l/rev vibratory 
response of the blades in the rotating system. 

In addition to the N/rev or NP vibration, there are secondary excitation frequencies 
which also must be taken into account. Any unbalance of the rotor system, such as blades 
out of track, will give rise to 1/rev or IP excitation of the fuselage. Dissymmetry of blade 
dampers may cause troublesome 2/rev excitation. Also important is the 1 x tail rotor 
rotational speed. 

Forward speed plays a critical role in vibrations. Figure 2.1 is a typical plot of the 
variation in vibration levels as a function of airspeed. The increase in vibrations levels in 
the 30 to 45 knot range can be attributed to blade vortex interaction due to the proximity 
of the rotor to its wake and the resulting non-uniform induced velocity field. 




Figure 2. 1 : Typical Helicopter Vibration Levels Versus Airspeed. After Ref. [2], 
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In the high speed flight regime, the increase in vibration with increase in airspeed is 
attributed to two sources. First, the aerodynamic excitation due to impingement of the 
rotor’s wake on the helicopter’s horizontal empennage. Second, amplification of blade 
section aerodynamic loads associated with increased cyclic pitch accompanied by 



increasing differences in dynamic pressure per blade revolution. Figure 2.2 shows the 
non-distorted rotor wake geometry for both low and high forward flight speeds. 




2. Random Aerodynamic Excitation 

This type of vibration is often associated with the turbulent downwash from the 
rotor, which impinges on the fuselage and its components. Similar to white noise, this 
excitation contains many harmonics that can excite fundamental fuselage or empennage 
modes. Often, an interesting characteristic of this type of excitation is that it is present 
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only within a certain range of airspeeds. Successful solutions of this problem have 
included modifying the structure excited and by altering the flow from the rotor. 

3. Self-Excited Vibrations 

A self-excited vibration is one that results in divergent oscillations where the 
system damping is negative. In the case of negative damping, the damping force, which is 
now a driving force, does positive work on the system. The work done by this force is 
converted into the additional kinetic energy of the increased vibration. A self-excited 
vibration cannot exist without an external source of energy. In the case of flutter, the 
airflow provides the source of energy. 

B. VIBRATION SPECIFICATIONS 

Vibration specifications for Army aircraft are governed by Aeronautical Design 
Standard (ADS) - 27 and MIL-STD-810D [Ref. 1], Human factors vibration 
requirements are expressed by ADS-27 in terms of an intrusion index which is the square 
root of a weighted sum of squares of four vibration components in each direction. The 
purpose of this is to bring together the twelve most significant crew vibration components 
into a single number for evaluation. A separate specification exists to capture the lower 
frequency one per main rotor revolution, or IP, vibration. 

The changes in vibrations that result for airspeed and maneuvering are addressed in 
ADS-27 by defining limits in terms of flight regions. There are three defined regions, with 
the first being' normal flight up to cruise speed, then maneuvers having a duration less than 
3 seconds, and lastly maneuvers having a duration greater than 3 seconds. There is also a 
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region for gunfire vibration. The specifications for the RAH-66 crew seat are listed in the 
Performance Weapon System Specification [Ref. 1] and are shown in Table 2.1. 
Additionally, Figures 2.3 and 2.4 display the currently measured intrusion index in the 
Pilot’s and Copilot’s crew station versus the specification. 





Intrusion Index 


MR IP (in/sec) 


Flight Region I 


1.2 


.15 


Flight Region II 


3.0 


.30 


Flight Region III 


4.0 


.40 


Weapons Firing Increment 


+2.0 


N/A 



Table 2. 1 : Crew Seat Intrusion Index and MR1P Vibration Limits. 



Because of the frequency bands that are used to calculate the intrusion index, any 
random, broad band, excitation occurring between main rotor IP and NP is not totally 
captured. Therefore, the intrusion index does not indicate the contribution of the broad 
band excitation to total response. 
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Figure 2.4: Intrusion Index - Copilot’s Station. 
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C. COMANCHE VIBRATION SOURCES 

1. Main Rotor Hub Forces and Moments 

As noted earlier, the primary source of helicopter fuselage vibrations is the rotor 
system. As the forces and moments are transmitted from the rotating to the fixed fuselage 
system, the rotor system tends to act as a filter. It can be shown for an “N” bladed rotor 
system that, assuming all blades see the same loading at the same azimuth position, only 
NQ frequency will be seen in the fixed system. This Nf2 frequency in the fixed system can 
arise from either N-l, N, or N+l excitation in the rotating system. Figure 2.5 is a 
schematic of the hub showing the mechanism of blade-fuselage coupling. 
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With the coordinate system given in Figure 2.5, forces in the rotating blade 
coordinate system can be related to forces in the fixed fuselage using the following 



transformation: 
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where [X,] is the transformation relating forces in the rotating r\, Q system to forces in the 
fixed X, Y, Z system. 

Thus for the five-bladed Comanche helicopter: 

• 4/rev and 6/rev flapwise blade root shears result in 5/rev hub pitching and 
rolling moments in the airframe. 

• 5/rev flapwise blade root shears feed into the fuselage directly as 5/rev vertical 
forces. 

• 4/rev and 6/rev chordwise root shears produce 5/rev airframe hub forces in the 
fore and aft and lateral directions. 

• 5/rev chordwise root shears results in 5/rev hub yawing moments. 
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For the RAH-66, the nominal rotor speed is 354.9 rpm. This places the 1/rev or 
IP excitation at 5.915 Hz, with the 5P blade passage frequency occurring at 29.575 Hz. 
[Ref. 3] 

2. Aerodynamic Excitation 

Aerodynamic excitation to the airframe and empennage arises from several 
sources: 

• Main rotor wake impingement. 

• Main and tail rotor downwash (induced flow) impingement. 

• Flow from around the main rotor hub and pylon assembly structure 

• Flow from attached aerodynamic surfaces. 

• Flow from protruding irregular fuselage structure, ranging from FABS 
(forward avionics bays) to as small as a canopy or door handles. 

The interaction of this flow is a complex subject, very difficult to predict, and the 
subject of numerous ongoing research projects. Rotorcraft can have problems where 
rotor airloads affect airframe loads and, in turn, the airframe effects rotor airloads and 
blade dynamic response. Some examples: The Sikorsky H-3 had a lateral tail shake that 
had its onset at about 120 knots. The problem was solved by installing a fairing or 
“beanie” on top of the main rotor hub to smooth out the flow from this region which was 
impinging on the aft fuselage. 

In development of the Boeing YUH-61 and Sikorsky YUH-60 helicopters, both 
companies obtained a significant improvement in N/rev vibration levels by “raising the 
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mast”, that is, increasing the separation between the main rotor hub and the fuselage. In 
development of the Hughes YAH-64, the initial T-tail design was scrapped and a movable 
horizontal stabilator installed to eliminate an excessive vibration problem associated with 
the T-tail. A visual summation of the typical sources of tail excitation is shown in Figure 
2 . 6 . 



Main Rotor Vortices 




Random Excitation 
from Flow Separation 



Figure 2.6: Typical Sources of Tail Excitation 

An experimental investigation of the interactions between a rotor and a T-tail 
empennage was conducted by Lieshman and Moedersheim [Ref. 4] at the University of 
Maryland Rotorcraft Center in 1997. Using a 1/6 scale four-bladed helicopter rotor with 
generic body shape, a T-tail empennage was placed down stream of the rotor hub as 
illustrated in Figure 2.7. Both the vertical and horizontal tail surfaces used the NACA 
0012 airfoil section with a constant chord of .203 meters. The horizontal tail was located 
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in two different positions, a high and low plane configuration, with zero degrees angle of 
attack. Additional parameters of the wind tunnel test are contained in Reference 4. 



Figure 2.7: Schematic of Rotor/Body/Empennage System. From [Ref. 4], 



Through their analysis, they drew the following conclusions: 

• Flow behind the rotor was highly asymmetric. The two low pressure regions 
trailing from the edges of the rotor disk provided evidence of the wake rolling- 
up into two larger vortex bundles. 

• At low advance ratios the flow over the horizontal tail was stalled. The flow 
gradually reattached as the advance ratio was increased. Airloading was 
asymmetric primarily due to the difference in strength of the two trailing vortex 
bundles trailing from the sides of the rotor disk. In some test conditions, a 
sharp trailing edge peak pressure was observed at the trailing edge of the 
horizontal tail. This indicated the potential of a junction vortex originating at 
the junction of the vertical and horizontal tail. 

• Measured time-varying pressures showed the flowfield near the horizontal tail 
was highly unsteady. Below advanced ratios of . 15 the unsteady pressure 
response was low. As the wake transitioned on to the tail, unsteady pressures 
increased significantly. 

• Unsteady pressure responses are characterized by two events, unsteady lift 
produced by the rotor wake and disturbances caused by the convection of 
individual vortices. 
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For the RAH-66, the rotor wake transitions above the T-tail between 40 and 60 
knots, as shown in Figure 2.5 [Ref. 3], Beyond this airspeed range the actual wake is not 
impinging on the horizontal stabilizer. From Leishman’s research; however, there can be 
convected vortices associated with the passage of blades that can induce disturbances and 
oscillating loading conditions. 




Figure 2.8: Rotor Wake Distance from Horizontal Stabilizer. After Ref. [3], 



In addition to the rotor vortices, the wake behind the rotor hub constitutes a low- 
pressure sink that can draw the flow off of the upper fuselage and produce separation 
[Ref. 5], The frequency content of this excitation is often determined by modeling the 
rotor hub as a circular cylinder where vortices are shed as a Karman vortex street. 
Through the Strouhal number, a non-dimensional constant, vortex shedding frequency can 
be related to ffeestream velocity and the diameter of a two dimensional cylinder. [Ref. 5] 
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The relationship is given by: 




su 

D 



where f s = Vortex shedding frequency (Hz) 

S = Strouhal number 
U = Freestream velocity (ft/ sec) 

D = Cylinder diameter (ft) 

Determination of the Strouhal number from Schlichting [Ref. 5:p:32] requires 
interpolating from a small plot. Using the results from a wind tunnel experiment done at 
the Naval Postgraduate School on a two-dimensional cylinder by Small, Hebbar, and 
Platzer [Ref. 6], the Strouhal number and vortex shedding frequency can more easily and 
accurately be determined. In this experiment [Ref. 6], the vortex shedding frequency and 
Strouhal number were calculated for Reynolds numbers ranging from 2.7x1 0 4 to 1.9xl0 5 . 
From this, the data presented in Figure 2.9 was generated. The Strouhal number 
represents the slope of the line and can be used to calculate the frequency of shedding 
from a know diameter and given velocity. 
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Velocity (ft/sec) 

Figure 2.9: Shedding frequency for the Quasi-two-dimensional Cylinder. After Ref. [6], 

For the Comanche, at 100 knots (168.78 ft/sec) and with a rotor hub pylon of 
approximately .75 meters (2.46 ft), the calculated vortex shedding frequency is 1 1.827 Hz. 
As the airspeed and point of separation changes, the shedding frequency will be effected, 
as shown by the equation on the previous page. The area around the main rotor pylon is 
actually a complex array of surfaces needed to reduce radar signature. Separated flow 
could occur from any of these surfaces, with varying frequency content of the vortices. In 
this case, more broadband vortex shedding will occur. 

From these factors, RAH-66 response data should show energy associated with 
main rotor 5P, 29.75 Hz and if a rotor track and balance problem, main rotor IP at 5.95 
Hz. There could also be energy associated with separated flow from the main rotor pylon, 
centered around 11-12 Hz, and some broad band excitation from the convected rotor 
vortices impinging on the horizontal tail. 
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III. SPECTRAL ANALYSIS 



Figure 3.1 shows the time history of lateral excitation at the pilot’s heels. When 
analyzed over the full record length of 60 seconds, the data is classified as random. Such 
random time functions as this data set commonly exist in nature; a typical example being 
pressure gusts aircraft encountered in flight. The method used to characterize these 
functions is known as spectral analysis, and a discussion of the theory, from Thomson 
[Ref. 7], and techniques for digital signal analysis follows. 



Lateral Response Pilot's Heels 

0.25 

1 °' 15 
O 0.05 

2 

jjj -0.05 
o 

o -0.15 

< 

-0.25 

0.0 0.5 1 0 1.5 2.0 2.5 3.0 3.5 4 0 

Time (sec) 

Figure 3.1: Time History of Lateral Acceleration. 




A. TIME AVERAGING 

Despite the irregular nature of the time history in Figure 3.1, most “random” 
functions exhibit some degree of statistical regularity and by applying averaging processes, 
we can characterize the phenomena. In observing the data, each record is referred to as a 
sample and the collection of samples is the ensemble. If the ensemble average is taken at a 
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specific time, ti , and then at another time, tj + r , and they remain the same, the process is 
said to be stationary. Similarly, if the ensemble averages are replaced by time averages, or 
averaging each sample, and the averages from different samples remains the same, the 
random process is ergodic. Key assumptions necessary for our analysis is that the random 
response phenomena exhibited by the aircraft is stationary and ergodic. As such, statistical 
properties can be derived from single time histories of sufficient length. [Ref. 7] 

To derive these properties a review of statistics is required. From the time 
histories, we can define the expected value (mean): 



For the case of discrete variables, the expected value is given by the following equation: 



We can then apply the averaging process to the mean square value, also referred to as the 
square of the expected value. 
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Linking this back to basic statistics, with the variance defined as: 




the mean square value (MSV) is simply the variance when the data has an expected value 
(mean) of zero. If thought of in terms of a sin wave, or multiple harmonic function, where 
the mean is removed (detrended) before analysis, then the variance and MSV are the same. 

Knowing these statistical properties of a given signal, the correlation, or similarity, 
between different signals can be determined. This is an important function in time domain 
analysis as it can help identify signals buried in noise or propagated in time through a 
structure. Additionally, it provides a transition into Fourier analysis and the frequency 
domain, where a structure’s spectral characters can better be analyzed. The correlation 
function is computed by multiplying the ordinates of two records at each time and 
determining the average value by dividing the sum by the number of products. The 
product is a maximum when the two signals are very similar or identical. For dissimilar 
records, the products will be negative and positive, so the sum will be small. If two 
signals are identical with one shifted by some time, we can compute the correlation by 
taking just a single time record and multiplying it by itself, with the specified time delay i. 
This is called the autocorrelation function, and is defined by the following formula. 
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If the time delay x = 0, the equation reduces to the mean square value or variance. 

R( 0) = x* = a 2 

If the two records are truly dissimilar, we call their product the cross correlation. It is 
define it by: 

R V ( r ) = £[*(0 y{t + T)] = ( X(t)y(t + r>) 

1 % 

= limf \xit)y{t + r)dt 

T ^oo d j/ 

/2 

B. FREQUENCY DOMAIN ANALYSIS 

Having defined averaging processes for the time domain, the random function can 
be transitioned to the frequency domain [Ref. 7], This will assist in relating the response 
data back to the sources of excitation discussed earlier. 

The time signal x(t) can be expressed in terms of its Fourier series: 

x(0 = fyy-' 

-00 

Because the record is periodic it can be represented by the Fourier series: 



x(0 = Re£O? 

n = 1 



2 
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where C„ a complex function and C„* is its complex conjugate. It can then be proven [Ref. 
7] that the mean square value of the function is given by 

?=lk,c; 

From this, it can be seen that the MSV can be computed in either the time or frequency 
domain, depending on which technique provides better analysis of the data. 

The mean square value (MSV) represents contributions from the entire frequency 
band of the signal. If only a certain frequency range is of interest, the MSV for that range 
of frequency, Af, is the power spectrum (PS). The power spectrum, P(f n ), is then just the 
product of the Fourier coefficients from the frequency range of interest,./*, with the mean 
square value being the sum of the power spectrum across the entire bandwidth. [Ref. 7] 

f(/.)=kc; F=i>(/„) 

^ n = 1 

If the power spectrum is divided by the frequency interval, Af, the result, defined 
as the discrete power spectral density(PSD) S(fJ, gives the contribution to the power 
spectrum from an even more narrow frequency band. This can be extended for a 
continuous spectrum as the limiting case of S(fJ as Af —> 0. 

P ff x Pjfn).C H C m 

Un) Af 2 A f 

The mean square value can then be computed from the PSD by 

n = 1 
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The root mean square (RMS), synonymous with the standard deviation, can be found by 
taking the square root of the mean square value (MSV). 



RMS = Jx T 



As we are dealing with non-periodic functions over discrete time intervals, use of 
the Fourier integral, a limiting case of the Fourier series, is required. The Fourier Integral 
and Fourier transform of x(t) are defined respectively as 



*(0 = f xtjy^df X(f) = j xtiy'^dt 

-oo -oo 

Using Parseval’s theorem, the MSV of the discrete signal in the time domain can 
be equated to the Fourier transform of the signal multiplied by its complex conjugate. If 
Xi(t) and x 2 (t) are real valued time functions and their Fourier transforms are given by 
X\(f) and X 2 (f) then 



J*. «*,«<*= \x,wcl<JW 



If X](t) = x 2 (t) - x(t) the MSV and PSD can be expressed in terms of the Fourier transform 



x 



2 



= lm7 /*’(')* =jlim7 

T-+co 1 _J/ T-+co 1 



-Vi 



X(J)X\f)df = \S(JW 



% 



24 



As with the Fourier coefficients, the MSV and RMS of the signal can be 
determined from the Fourier transforms. Additionally, the autocorrelation function, 
symmetric about the origin, can be related to the spectral density function via the Weiner- 
Khintchine equation. [Ref. 7] 

oo 

S(f) = 2 ^R{r)cos, 27 rfrdr 

o 

The cross-spectral density (CSD), defined as the Fourier transform of the cross- 
correlation function: 

V/)= ]R,Xr)e-™dT 

-00 

Can also be determined by multiplying the Fourier transforms of the two signals of 
interest. 

V/> = lim 7*W(/) 

T — >oo 1 

Given all of these tools, the one of major interest in analyzing the flight test data is 
the power spectral density (PSD). This function of frequency gives an indication of the 
contribution to the total response, in our case accelerations, due to a discrete frequency 
bandwidth. Additionally, the cross-spectral density recaptures the phasing information 
lost when the Fourier Integral is multiplied by its complex conjugate. 
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c. 



DIGITAL SIGNAL PROCESSING 



Given this theory, digital signal processing methods have been devised that 
compute the PSD from discrete time histories. Vectors represent the time histories. In 
our case, the sampling rate is 300 samples per second and record length is 60 seconds. 
Therefore, the vector of the discrete signal has 18,000 data points. 

To convert this time signal into either a PSD or CSD, the first step is to decide 
how many points, N, are needed in the Fast Fourier Transform (FFT). For computational 
speed and accuracy, FFTs are done in power of two. The number of points also 
determines what the effective bandwidth of the processed information will be. Bandwidth, 
Af, is defined as 1/Time, and for a discrete signal, Time = N x At (time increment). The 
time increment, At, is equal to 1/sampling rate, or, in our case , 1/300 sec. If a 1024 
point FFT is utilized, Af = 300/1024 = .293 Hz. As the number of points in the FFT 
increases, the bandwidth decreases and you get finer frequency resolution. 

There is a penalty when working with random data to increasing the number of 
points in the FFT. Based on N, the signal is divided into rid records or windows and has 
any linear trend removed. As N increases, the number of windows, n<j, decreases and there 
is less averaging occurring. This is demonstrated when computing the expected value. 



% 
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To suppress leakage of the signal from one frequency range to the next, time windows are 
used that taper the time history data. This eliminates discontinuities at the beginning and 
end of the record to be analyzed. A commonly used window, that is used exclusively 
during this analysis, is the Hanning window. It is given by 

2 izt 

w(t) = .5 * (1 - cos(— )) 

and multiplied by the N data points of each nj segment. Use of the window does reduce 
the amplitude of the processed signal. In the case of the Hanning window, the output is 
.375 times smaller. A correction is made in the normalizing of the PSD to account for 
this. 

Variance is introduced by tapering the signal. To reduce this, the n<j segments are 
overlapped. In terms of the Hanning window, the overlapping brings the "peaks" of each 
cosine window closer to each other. This process, basically, takes the points that were at 
the end of one n<j segment, where the window is near zero, and makes it the "main area" 
for the next. The FFTs are then computed for each tapered, overlapped segment. 

From the FFTs, the PSD can be determined using the following formulas: 

±(\fft\X 2* E[<\FFT\ 2 )]* M 

^ “ ( normiwindow )) 2 

The expected value is the sum of the FFTs of each segment multiplied by its 
complex conjugate, or the magnitude of the FFT squared. Either way, it produces a real 
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valued function that is then divided by the number of segments, ru As discussed earlier, 
and seen here, increasing N reduces the number of FFT averages in the expected value. 

The PSD is then determined by "normalizing" the expected value. This first 
involves multiplying by 2 as the PSD represents the one-sided autospectral density 
function. Then, the signal is multiplied by the time increment, At (or dividing by the 
sampling rate, from the definition of At) and divided by the norm of the window, squared. 
The reasons for this come from the definition of the FFT and by canceling several terms. 
Bendat and Piersol [Ref. 8:p.371] provide the exact proof; however, their notation 
requires dividing the signal by N and then adjusting for the attenuation caused by a 
tapering window. The norm of the window (its length in N dimensions) accounts for both 
elements. For no windowing, computationally done using a boxcar window, the norm 2 is 
the number of points in the FFT, 1024. For a Hanning window, the norm 2 = 1024 x .375 
= 384. As discussed earlier, this corresponds to the attenuation causing by using a 
Hanning window. Figure 3.2 and 3.3 provide a visual illustration of the procedure used to 
compute the autospectral density function (PSD). 

For the cross spectral density, the same steps are followed for each time history. 
Then the FFTs from each signal, x(t) and y(t) are multiplied by each other and the 
expected value is a complex function. E[FFT *] is normalized using the same method. 
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The time history is divided into segments. For computational speed and accuracy, each segment 
consists of 2 n data points. The data is detrended to eliminate the stationary component. 
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A tapenng window, such as a Hanning, can be used to reduce the spectral leakage. To 
counter the variance caused by tapering, overlapping time histories are used. 
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Figure 3.2: PSD Computational Steps - Time Domain. 
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The FFT is computed of each tapered time history. The magnitude of the FFT is then squared. 
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Figure 3.3: PSD Computation - Frequency Domain. 
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The PSD or CSD output vector is of length N/2 with a frequency band from 
0:Af:(sampling rate/2). As an example, with N of 1024 points and a sampling rate of 300 
Hz, the maximum frequency, commonly referred to as the Nyquist frequency, is 1 50 Hz. 
The frequency vector begins at 0 and goes to 150 Hz in increments of 300/1024 or .293. 
For each point of the frequency vector, the PSD vector has a corresponding value. This 
value represents the energy associated with that frequency band. So, in this example, the 
first value of the PSD vector represents the energy present within the 0 to .293 Hz range. 

Units of either the PSD or CSD are (input units) 2 /Hz. In the case of accelerometer 
output, the PSD units are g 2 /Hz. With units like this, it is often difficult to grasp an 
understanding of what the information really means. The RMS is much easier to interpret, 
which has the same units of the input, in this case, g’s. To compute the RMS, the PSD 
vector is multiplied by the frequency increment, Af. This vector is then summed within the 
frequency range of interest. The resulting scalar is the MSV and the RMS is the square 
root. If, for example, the RMS between 5-25 Hz is needed, the PSD vector is multiplied 
by Af = 300/1024. Then, only the indices of the vector that represent energy between 5- 
25 Hz are summed, points 17-86 of the 512 total PSD values. The square root is the 
RMS and represents the total energy or response between 5 and 25 Hz. Comparing back 
to the time domain, this RMS value equates to the standard deviation. 
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IV. FLIGHT TEST DATA ANALYSIS 



Time histories of accelerometers and pressure transducers from flight tests were 
analyzed using both 1024 and 4096 point FFT. For computational speed, 1024 point 
FFTs were utilized for final analysis. Several different tapering windows were also applied 
with no noticeable difference. The excitation frequencies are clearly evident and are an 
order of magnitude higher than the energy of the surrounding frequency bands. Therefore, 
even when no side lobe reduction techniques were utilized, the frequency content of the 
response was pronounced. To remain consistent with the methods used by Sikorsky, a 
Hanning window with 67% overlap was used. This did reduce the spectral leakage and 
variance associated with windowing and produced a more accurate PSD value. 

A. ACCELEROMETERS 

Data records at 60, 80, 100, 120, and 140 knots were obtained for the 
accelerometers listed in Table 4.1. The sampling rate for all signals was 300 samples per 
second with a 60 second record length. A complete listing of all test flight parameter is 
contained in Appendix A. 

The PSD of each accelerometer through the range of airspeeds was computed. In 
looking at the energy across the entire frequency spectrum, 0-150 Hz, there was 
significant energy from 0-30 Hz, and none of magnitude beyond that range. Figure 4.1 
through 4.3 display the PSD for nose cone lateral, pilot heel lateral, and copilot heel 
vertical. The programs used to compute the PSDs and CSDs is in Appendix B. 
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Accelerometers 


Mnemonic 


Measurement 


LNOSE:D:300 


Lateral Nose 


VNOSE:D:300 


Vertical Nose 


LONOSE:D:300 


Longitudinal Nose 


LPHEEL:D:300 


Lateral Pilot Heel 


VPHEEL:D:300 


Vertical Pilot Heel 


LPF:D:300 


Lateral Pilot Floor 


VPF:D:300 


Vertical Pilot Floor 


LOPF:D:300 


Longitudinal Pilot Floor 


LPOVRHD:D:300 


Lateral Pilot Overhead 


VCOPHEEL:D:300 


Vertical Copilot Heel 


LCOPHEEL:D:300 


Longitudinal Copilot Heel 


VCOPF:D:300 


Vertical Copilot Floor 


LCPOVRHD:D:300 


Lateral Copilot Overhead 


LVPYLNT:D:300 


Lateral - Vertical Pylon Top 


VSTABTLT :D:300 


Vertical - Hor Stab Left Tip 


LOSTABTL:D:300 


Longitudinal - Hor Stab Left Tip 


VSTABTRT:D:300 


Vertical - Hor Stab Right Tip 


LOSTABTR:D:300 


Longitudinal - Hor Stab Right Tip 



Table 4. 1 : Accelerometers. 



As previously noted, primary excitation for a 5-bladed rotor occurs at 5P, which is 
29.6 Hz for the Comanche. As expected for this case, 5P response centered at 29.6 Hz 
dominates, especially at higher airspeeds. In contrast, the amplitude of the other frequency 
bands is not discemable. To provide more insight at response due to other than 5P, 
Figures 4.4 through 4.18 displays the accelerometer response from 5-25 Hz. Within this 
frequency band, much greater detail in the response data can be discerned. Given on the 

next 1 7 pages are plots that indicate the frequencies of excitation along with the expected 

* 

increase with airspeed. 
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g /Hz 



PSD Nose Cone - Lateral Response 




Airspeed Frequency 

Figure 4. 1 : Nose Cone - Lateral 5-30 Hertz. 
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g /Hz 



PSD Pilot Heel - Lateral Response 




Airspeed Frequency 



Figure 4.2; Pilot Heel - Lateral 5-30 Hertz. 
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g /Hz 



PSD Copilot Heel - Vertical Response 




Airspeed Frequency 

Figure 4.3: Copilot Heel - Vertical 5-30 Hertz. 
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^ q -3 PSD Nose Cone - Vertical Response 
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Figure 4.4: Nose Cone - Vertical 5-25 Hertz. 
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PSD Nose Cone - Lateral Response 
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Figure 4.5: Nose Cone - Lateral 5-25 Hertz. 
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^q- 3 PSD Pilot Heel - Lateral Response 




Airspeed Frequency 



Figure 4.6: Pilot Heel - Lateral 5-25 Hertz. 
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g /Hz 



1 q- 3 PSD Pilot Heel - Vertical Response 




Airspeed Frequency 



Figure 4.7: Pilot Heel - Vertical 5-25 Hertz. 
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^ 0 3 PSD Pilot Floor - Lateral Response 




Airspeed Frequency 



Figure 4.8: Pilot Floor - Lateral 5-25 Hertz. 
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1 q -4 PSD Pilot Floor - Vertical Response 




Airspeed Frequency 

Figure 4.9: Pilot Floor - Vertical 5-25 Hertz. 
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ig-3 PSD Copilot Heel - Vertical Response 
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Figure 4.10: Copilot Heel - Vertical 5-25 Hertz. 
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Airspeed Frequency 

Figure 4. 1 1 : Copilot Floor - Vertical 5-25 Hertz. 
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PSD Pilot Overhead - Lateral Response 




Airspeed Frequency 



Figure 4.12: Pilot Overhead - Lateral 5-25 Hertz. 
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PSD Copilot Overhead - Lateral Response 
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Figure 4.13: Copilot Overhead 



Frequency 
Lateral 5-25 Hertz. 
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PSD Vertical Pylon Top - Lateral Response 




Airspeed Frequency 



Figure 4. 14: Vertical Pylon Top - Lateral 5-25 Hertz. 
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g /Hz 



PSD Horizontal Stab Left Tip - Vertical Response 




Airspeed Frequency 

Figure 4.15: Horizontal Stab Left Tip - Vertical 5-25 Hertz. 
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PSD Horizontal Stab Right Tip - Vertical Response 




Airspeed Frequency 



Figure 4.16: Horizontal Stab Right Tip - Vertical 5-25 Hertz. 
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g /Hz 



PSD Horizontal Stab Left Tip - Longitudinal Response 




Airspeed Frequency 

Figure 4.17: Horizontal Stab Left Tip - Longitudinal 5-25 Hertz. 
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g /Hz 



PSD Horizontal Stab Right Tip - Longitudinal Response 
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Figure 4.18: Horizontal Stab Right Tip - Longitudinal 5-25 Hertz. 
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From a review of the plots, there are three frequency bands that continuously 
display energy. Energy in this context is based upon measured vibration response data as 
either acceleration in g’s or PSD in (g’s) 2 /Hz. The first band is centered around 5-8 Hz, 
the second from 10-13 Hz, and the last 22-25 Hz. From the PSD data, the RMS for each 
airspeed was calculated. Figures 4.19 through 4.22 show the computed RMS from the 
nose cone, pilot’s station, and copilot’s station. The RMS including and not including 5P 
is plotted along with both lateral and vertical response. The 5P vibration adds a significant 
amount to the total response and produces a sharp rise in accelerometer readings 
acceleration at airspeeds above 100 knots. The total response for the frequency band of 
5-25 Hz begins to level off at 100 knots, the slope being much less than the RMS with 5P 
included. This excitation, not including 5P, is less prevalent at lower airspeeds and 
significantly increases from 80-100 knots. Past 100 knots, this frequency band energy 
begins to level off while 5P vibration levels increase significantly. By plotting both lateral 
and vertical RMS, comparisons as to the more dominant response can be made. Lateral to 
vertical ratios typically range of the order from 1.25:1 to 2.0:1. In all cases except for the 
pilot’s floor, lateral response dominates. To capture the maximum excitation outside of 
the 5P range, and minimize 5P response, a more detailed analysis is presented of the PSD 
at 100 knots. 

Figures 4.23 through 4.28 present the PSD for 100 knots. Again, both lateral and 
vertical responses are plotted together. As was seen in the RMS calculations, the lateral 
response is dominant. Additionally, the three frequency regions identified earlier are 
evident. 



53 



Acceleration (g's) Acceleration (g s) 



Root Mean Square Nose Cone 




Figure 4. 19: Root Mean Square - Nose Cone. 



Root Mean Square Pilot Heels 




Figure 4.20: Root Mean Square - Pilot Heels. 
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Acceleration (g s) Acceleration (g's) 



Root Mean Square Pilot Floor 




Figure 4.21: Root Mean Square: Pilot Floor. 
Root Mean Square Copilot Station 




Figure 4.22: Root Mean Square - Copilot Station. 
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Frequency (Hertz) 

Figure 4.23: Nose Cone - 100 Knots. 
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Figure 4.24: Pilot Heels — 100 Knots. 
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PSD Pilot Floor - 1 024 Point FFT 
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Figure 4.25: Pilot Floor - 100 Knots. 
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Figure 4.26: Copilot Station - 100 Knots. 
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g 2 /Hz g /Hz 



PSD Horizontal Stab Lt Tip - 1024 Point FFT 
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Figure 4.27: Horizontal Stab Left Tip - 100 Knots. 
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Figure 4.28: Horizontal Stab Right Tip - 100 Knots. 



58 



To quantify the amount of response due to these three frequency bands, or bins, 
RMS values for the discrete ranges of 5-8 Hz, 10-13 Hz, and 22-25 Hz were calculated. 
Each RMS value was summed and compared to the total RMS, not including 5P response. 
The data is presented in Table 4.2. Although these three frequency bins represent only 9 
Hz, or 45% of the total bandwidth, they constitute from 50% of the total response to as 
much as 90%. Comparing the bins, the more dominant frequency range was 10-13 Hz, 
with 60% of the accelerometers having their peak response in this region. 



Root Mean Square - 1 00 Knots - q's 




Frequency Bins 


% Total 
RMS 


Accelerometer 


5-8 Hz 


10-13 Hz 


22-25 Hz 


Nose Cone Vertical 


0.0214 


0,0310 


0.0288 


51.87 


Nose Cone Lateral 


0.0448 


0.0920 


0.0306 


91.54 


Pilot Heel Lateral 


0.0352 


0.0573 


0.0159 


92.54 


Pilot Heel Vertical 


0.0161 


0.0283 


0.0096 


57.09 


Pilot Floor Lateral 


0.0238 


0.0191 


0.0058 


92.73 


Pilot Floor Vertical 


0.0095 


0,0155 


■ 


0.0050 


64.88 


CPG Heel Vertical 


0.0054 


0.0103 


0.0147 


82.59 


CPG Floor Vertical 


0.0030 


0.0052 


0.0162 


65.97 


Pilot Overhead Lat 


0.0215 


0.0816 


0.0868 


88.90 


CPG Overhead Lat 


0.0061 


0.0229 


0.0647 


02 ♦$4 | 


Vert Pylon Top Lat 


0.4521 


0.1669 


0.1024 


87.02 


HorStab Lt Tip Vert 


0.2928 


0.5035 


Q„33§5 


84.60 


HorStab Rt Tip Vert 


0.4284 


0.7390 


0.5327 


81.18 


HorStab Lt Tip Long 


0.2325 


1.0588 


0.3360 


75.25 


HorStab Rt Tip Long 


0.2354 


1.0268 


0.0987 


86.69 



Table 4.2: RMS of Frequency Bins - 100 Knots. 
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B. FLIGHT TEST PRESSURE MEASUREMENTS 



The differential pressure transducer data records that correspond to the 
accelerometer data are listed in Table 4.3. The sampling rate remained 300 samples per 
second with a 60 second record length. 



Differential Pressures 


Mnemonic 


Location 


VERTDP2 


Differential Pressure, 7% chord, 30% 
up Vertical Tail 


VERTDP4 


Differential Pressure, 7% chord, 60% 
up Vertical Tail 


LHORZDP3 


Differential Pressure, 7% chord, 60% 
Span Left Horizontal Tail 


RHORZDP4 


Differential Pressure, 7% chord, 60% 
Span Right Horizontal Tail 


Table 4.1 


1: Differential Pressure Transducers. 



The frequency content of the differential pressures was determined and compared 
to the response PSDs, findings of Leishman and Moedersheim [Ref. 4], and the 
approximated shedding frequency of the rotor hub pylon. The analysis also provided 
insight into which areas of the structure are excited. This aided in focussing the 
NASTRAN analysis, deriving the transfer functions between the random aerodynamic 
loading and the response locations, and correlation with mode shapes. This will be 
addressed in Chapter VI. Figure 4.29 shows the variation of differential pressure from the 
left horizontal stabilizer as a function of airspeed from 60 to 140 knots. 



60 



PSD Diff Pressure 7% Chord 60% Span Lt Horizontal Stab 




Airspeed Frequency 

Figure 4.29: Differential Pressure - Left Horizontal Stabilizer. 



As with the response data, the magnitude of the 5P differential pressure (29.85 Hz) 
is an order of magnitudes greater than excitation in any other frequency band. To gain a 
better appreciation for the vibratory pressures at other than 5P, the transducer output is 
plotted for the 5-25 Hz range in Figures 4.30 through 4.33. 

Given are three dimensional plots of airspeed, frequency and magnitude of the 
PSD. The asymmetry between the left and right horizontal tail is expected due to the 
angle of the vertical tail. The right horizontal tail exhibited pressure spikes at the 
previously noted discrete frequency bands of 5-8 Hz, 10-13 Hz, and 22-25 Hz. 
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Figure 4.30: Differential Pressure - Right Horizontal Tail. 
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Figure 4.31: Differential Pressure - Left Horizontal Tail. 
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Figure 4.32: Differential Pressure - Vertical Tail. 
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Figure 4.33: Differential Pressure - Vertical Tail. 
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Unlike the discrete excitation of the right horizontal tail, the left horizontal tail has 
broad band excitation from 11-19 Hz. Both span locations on the vertical tail show 
symmetry in the shape of broad band excitation, similar to that of the left horizontal. 

There is a decrease in magnitude, on the order of 12: 1, as you move up the vertical tail, 
from the 30% to 60% span location. 

Figures 4.34 through 4.35 are the PSDs of differential pressures at 100 knots. 

Both horizontal and vertical transducers are contained in the same plot for better 
comparison. For the horizontal tail, the magnitude of the pressures seen on the left is 
much great than that of the right. The energy is peaked within the 1 1 .5 to 14 Hz band, 
with an additional peak at 18 Hz. The right displays its peak pressure at 6 Hz, which 
corresponds with the IP of the aircraft. 

The vertical tail, at both the 30% and 60% span location, have the same 
frequency distribution as the left horizontal, 1 1 .5 to 14 Hz. As seen earlier, the magnitude 
of the pressure continues to diminish as you move up the vertical tail. Although the 
pressure distribution between the vertical tail and left horizontal is similar, the pressures 
seen at 30% span of the vertical are 3.5 times greater than that of the left horizontal. The 
pressures at the 60% span of the vertical and the left horizontal are of equal magnitude. 
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Figure 4.34: Differential Pressure Horizontal Tail - 100 Knots. 
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Figure 4.35: Differential Pressure Vertical Tail - 100 Knots. 
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Linking the pressure distributions back to response data, RMS values were 
calculated for the discrete frequency bins used earlier. The unit for the RMS of the 
pressures is psi. Table 4.4 presents both the magnitudes of pressures as you move up the 
vertical tail and along the horizontal as well as the dominant frequency band. 



Root Mean Square - 100 Knots - psi 




Frequency Bins 


Pressure Transducer 


5-8 Hz 


10-13 Hz 


22-25 Hz 


7% Ch 30% Sp Vert 


0.0263 


0.0503 


0.0268 


7% Ch 60% Sp Vert 


0.0182 


0.0206 


0.0184 


7% Ch 60 % Sp Lt Hor 


0.0121 


0.0268 


0.0135 


7% Ch 60% Sp Rt Hor 


0.0066 


0.0057 


0.0048 



Table 4.4: RMS of Differential Pressures - 100 Knots. 



Understanding the flow over the empennage, even for the limited amount of 
transducers, assisted in focussing the analysis of the structural response. The following 
chapter presents the structural analysis and correlation between aerodynamic forces acting 
on the surface and response at accelerometers. 
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V. STRUCTURAL DYNAMICS 



The dynamics of the multiple degree of freedom (MDOF) Comanche structure 
govern the transfer of random loads to fuselage response. In order to understand the 
MDOF system, the single degree of freedom case, also called a single modal oscillator, 
must first be studied. From these solution methods, the dynamics of the complex structure 
can then be quantified using a systematic modal analysis approach. Of particular interest is 
the response of the system to harmonic, or frequency dependent excitation. The quantity 
that describes this is called the frequency response function (FRF) or, from the Laplace 
domain, the transfer function evaluated at (/'ey). Working from a vibrations perspective, 
the FRF for a single DOF system will be computed and then applied to a MDOF 
Comanche structure. 

A. SINGLE DEGREE OF FREEDOM (SDOF) OSCILLATORS 

The response of a single degree of freedom system to a forced harmonic excitation 
can be determined from the equation of motion that describes the system. The SDOF 
oscillator is characterized by 

mx + cx + /ac = F cos(Q/) 

where 

m Total (modal) mass participating in mode of vibration 

k Total (modal) stiffness participating in mode of vibration 

c Total (modal) damping participating in mode of vibration 
F Magnitude of forcing function 
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The equation can then be divided by m and is in the generalized form. 



.. c . k F 

X H XH X = — COS(Q/) or 

m mm 



, F 

X + 2^(0 n x + co n x- — cos (Clt) 

m 



The total solution is comprised of two parts: 

*(0 = *homo(0+Vt (0 

The homogeneous solution, also known as the transient or free response is 

*homo (0 = e~^(C x cos (co d t) + C 2 sin(^/)) 
where the constants C/ and C 2 are determined by applying initial conditions. The 
particular solution is found for a specified forcing function. In the case of a cosine forcing 
function, the solution is 

x part (0 = ^cos(n/-O) 



and the system responds at the frequency of excitation, Q. 

To solve for the amplitude and phase of the response, A and (p, complex algebra is 
used and the general forcing function becomes 

F(t ) = Fe lCu = F(cos(Qt) + j sin(0/)) where F = — 

m 

With this definition, the general equation for the forced response becomes 

x + 2^a> n x + ar n x = Fe jQl 

-% 

and if the assumed solution for the particular solution is 

*pan (0 = AejQl 
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the necessary derivatives can be taken to find A in terms of known quantities 



(-Q 2 +2<£lcoJ + co 2 n )A = F 
2 

dividing by the natural frequency squared, o n and grouping the real and imaginary parts, 

A = Z(Q)*A = F 



f a 2 






+ 2C—J 


LI 


J 



with the force amplitude redefined as 

F F =■ 

-— = — => F 
k 

The terms inside the brackets represent the complex frequency dependent constant, 
Z(Q), known as the impedence, which is a function of the frequency fl The equation can 
also be written in the form 



Converting to polar form and dividing magnitudes and subtracting angles, a solution for A 

[z,^n)+jZ^(n)]A = F 



is obtained 



A = 



F 
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The modal parameters defined in our equation of motion can be substituted in to yield 



A = ■ 



F 



re ^ where <3> = tan' 1 
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The particular solution is then defined as 

* P an (0 = \A\e j(n, ~^ = |^|((cos(n? -<&)+j sin (Qi - O)) 

As the forcing function was a cosine function, only the real part is taken and the 
general solution is 

Xpat (0 = H(C0S(O/ - O) 



The amplitude of the steady state response is given by 




F 



1- 



n 



2 V 



CD 



n J 



2C — 



V 



CD, 



n J 



where 



F 



F 

k 



If the amplitude is then normalized by the static displacement F/k, the frequency response 
function (FRF) is found 



H( Q) = 



1- 



n 



2 A 



CD 



n J 






Q 



CD 



nj 
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The FRF can then be used to relate any forcing function at a given frequency to the 
output. If the forcing function is defined as F(fY) and the output, say displacement, is 
Y(F2), then the relationship between the two is given by 

7(Q) = H(Q)F(D.) 

B. MODAL ANALYSIS 

Knowing how to solve for the FRF of a SDOF system, these methods can be 
applied to the MDOF systems used to represent complex structures. In order to solve the 
coupled set of equations associated with MDOF systems, they first must be decoupled into 
single DOF modal oscillators. Each modal oscillator is then solved in the procedure just 
described and then combined to form the total solution. This method is called mode 
superposition and is consistent with the matrix methods currently used for most structural 
analysis. [Ref 9] 

If the MDOF structure is defined as a set of NxN coupled equations represented 
by 

[M]{x} + [C]{i} + [KKx}={F(,)} 

The solution to this equation is first found by solving the homogeneous eigenvalue 
problem with no damping, C. Substituting an assumed solution of the form {x}={4>}Oe' air 
and dividing by the mass, the eigenvalue problem is defined as 
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The non-trivial solution of the eigenvalue problem produces a vector of length N 
that contains the natural frequencies squared, ® n 2 with units of (rad/sec) 2 . The mode 
shapes, or eigenvectors {<)>}, can be easily solved for and combined into the modal matrix. 

[®] = [{®' Jj® 1 }{® s }..J® N H 

The modal matrix, although not an orthogonal matrix, is orthogonal with respect 
to the mass and stiffness matrix. As such, the mass and stiffness matrices can be 
transformed by the modal matrix to diagonal matrices that decouple the equations of 
motion and yield N - SDOF modal oscillators. The coordinate system must also be 
transformed to modal coordinates, q 

M=MM 

Pre-multiplying by the transpose of the modal matrix and substituting for x, the equation 
of motion, minus the damping becomes 

[®hmi®m+ [* nKMqH®r M 

Defining the modal mass, stiffness, and force as 

[®] T [mI®] = [m]; [®] T [kJ®] = [k]; {f}= [®f{F} 
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The modal equation of motion is 

[Mfc}+[£fe={p! 

These diagonal matrices now represent the N modal oscillator SDOF equations of motion. 
A solution for each can easily be obtained. 

Damping is re-introduced on a modal basis. As each differential equation is 
solved, the damping associated with that particular natural frequency, or mode, is applied. 
This is consistent with the techniques used in parameter identification where damping 
ratios are determined on a modal basis. The response in physical coordinates is found by 
multiply each modal solution by the modal matrix 



If interested in the frequency response function (FRF) between only certain points, 
i and j, on an MDOF system, the following method is used 



"p" modes 

H„(Q)= £ 

r = l 



O r O r 

> i 



<o r -0 2 +j2C r co r Cl 
Once the FRF is known for those points, the magnitude of the response, A, due to 

a harmonic excitation of amplitude F and at a driving frequency of Q, can be determined 

by multiplying the FRF and F. 
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C. RIGID BODY MODES 



Unconstrained systems have rigid body modes. This occurs when the structure 
translates or rotates without any deformation. Numerically, this appears in the solution as 
natural frequencies, co n , of zero. The reason for this is based on the relationship between 
the stiffness matrix, K, and the strain energy of the system. If the MDOF system moves 
without deforming, no strain energy is produced. The stiffness matrix, K, becomes 
positive semi-definite with a determinate of zero. As such, the solution to the eignevalue 
problem has as many trivial solutions, zero, as there are rigid body modes. For an aircraft, 
there are six rigid body modes, which is the case for the dynamics of any body in three- 
dimensional space. 

D. FINITE ELEMENT MODELING 

To determine the frequency response functions using modal analysis, the 
Comanche structure was modeled using the software package NASTRAN. The National 
Aeronautics and Space Administration (NASA) initially development NASTRAN in the 
1960s. NASTRAN is an acronym for NASA STRuctural ANalysis. Originally written in 
FORTRAN, it uses the finite element method (FEM) to analyze structural models [Ref. 
10], The foundation of this is the discretization of the structures stiffness and mass. The 
FEM provides the basis for algorithms that can efficiently analyze complex structures such 
as the Comanche. A detailed discussion of the Finite Element method and its relationship 
to structural dynamics is contained in Reference 1 1 . NASTRAN version 70 was used for 
this analysis. 
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The results generated using NASTRAN were viewed using PATRAN 7.0, an 
integrated computer-aided engineering (CAE) tool. The menu-driven graphical interface 
makes model manipulation easy and aids in visualization of mode shapes and model 
geometry. 
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VI. RAH-66 NASTRAN ANALYSIS 



A NASTRAN finite element model of the RAH-66 was provided by Sikorsky. This 
model, yrahj 1 57_reinf_2kg.dat, was modified by engineers at Boeing to replicate the 
configuration of the prototype aircraft. The name of the model is very important as 
Sikorsky manages many configurations and models, each with different stiffness and mass 
properties. The model used in this analysis included additional stiffening to the tail of the 
aircraft and 2.2 kg tip weights on both sides of the horizontal tail. With this, the natural 
frequencies, mode shapes, and frequency response functions (FRF) could be determined 
that most closely represented the actual aircraft. This was critical as these FRFs were later 
used with accelerometer and pressure transducer data to compute estimates of the loads. 

A PATRAN representation of the NASTRAN model is shown in Figure 6. 1 . 

A. NATURAL FREQUENCIES 

Hover frequency response testing of the test flight aircraft was conducted at West 
Palm Beach to support aeroservoelastic analysis. The purpose was to identify actual 
natural frequencies and modal damping. Through flight control frequency sweeps from 5- 
15 Hz, actual fuselage natural frequencies and percent critical damping were determined 
for selected modes. 
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Rather than adjusting the stiffness and mass properties of the NASTRAN model to 
return the actual natural frequencies, engineers at Boeing demonstrated that if the mode 
shapes from the original NASTRAN model were used along with adjusted natural 
frequencies, the resulting displacements and accelerations matched actual test flight 
response data. The method used to accomplish that was to first run the NASTRAN 
analysis deck. Then, during a restart using DMAP commands, the eigenvalue matrix 
(representing NASTRAN natural frequencies) was altered to actual frequencies. The 
resulting accelerations from application of harmonic loads closely matched actual response 
due to harmonic sweeps of the tail fan. Table 6. 1 displays the NASTRAN natural 
frequencies and actual natural frequencies and damping, for selected modes, as determined 
from flight test. 
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, . . ; 

Mode 


NASTRAN 
Freq (Hz) 


. . , . | 

Description 




tTest 
Damping 
% Critical 


7 


6.32 


Aft Fuselage Lat/Torsion 


6.00 




8 


7.77 


Fuselage Vertical Bending 


7.70 


2.34% 


9 


8.53 


Pilot Vert/Pitch 






10 


8.70 


Pilot Lat/Yaw 






11 


10.24 


Horiz Stab Roll 


9.50 


2.90% 


12 


10.88 


Copilot Vert/Pitch 






13 


10.98 


Copilot Lat/Y aw 






14 


12.84 


Horiz Stab Yaw 


11.60 


3.33% 


15 


13.67 


Tail Rotor Drive Shaft Mode 






16 


14.52 


2 nd Fuselage Vertical 






17 


14.99 


Quill Shaft 






18 


15.75 


Tail Rotor Drive Shaft 






19 


15.95 


Fuselage Head Vertical 






20 


16.07 


2 nd Fuselage Lateral 






21 


16.55 


Pilot Seat Lateral Roll 






22 


16.84 


Main Xmsn Lat/Roll 






23 


17.10 


Tail Rotor Drive Shaft 






24 


17.82 


Tail Rotor Drive Shaft 






25 


18.24 


Ammo Bay 






26 


18.82 


Copilot Lat/Roll 






27 


18.89 


Main Xmsn Lat/Roll 






28 


20.45 


Fuselage Mode 






29 


21.12 


Pilot Seat 






30 


21.56 


Stab Vertical Bending 






31 


21.78 


Tail Rotor Driveshaft 






32 


21.80 


Head Lateral 






33 


21.87 


Tail Rotor Driveshaft 






34 


23.05 


Copilot Sidearm 






35 


24.02 


Weapons Door 






36 


24.04 


Fuselage 






37 


24.23 


Weapons Door 






38 


24.53 


ECU Mode 






39 


24.61 


Pilot Seat/Sidearm 






40 


24.74 


Canopy 







Table 6. 1 : Aircraft Natural Frequencies. 
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B. MODE SHAPES 



Selected mode shapes for free response from the original NASTRAN solution are 
contained in Appendix C. Again, as the flight test derived natural frequencies were only 
slightly different from the NASTRAN frequencies, it is assumed that the NASTRAN 
mode shapes are accurate. Three critical mode shapes are shown in Figure 6.2 through 
6.5. The first two mode shapes were identified in test flight and have large lateral 
response in frequency bins previously discussed. The third mode shape is outside of the 
frequency range; however it correlates well with the response seen in the 22-25 Hz range. 

Figure 6.2 is the aft fuselage lateral torsion mode. This motion is consistent with 
the accelerometer output in the 5-8 Hz bandwidth. Referring to Table 4.2, correlation can 
be made between the flight test output and the free response. There is large response in 
this frequency bin from the accelerometer at the top of the vertical stabilizer. This same 
motion occurs in the 6 Hz mode shape. 

The horizontal stabilizer yaw mode at 1 1.6 Hz, displayed in Figure 6.3, has 
significant lateral response at the base of the vertical tail and at both the pilot station and 
nose cone. The RMS of accelerometer output from both the pilot station and nose cone 
have the largest magnitude in the 10-13 Hz bandwidth. 

There is no mode shape within the 22-25 Hz bandwidth that correlates well with 
the flight test data. Figure 6.4 shows the horizontal stabilizer bending mode at 26.02 Hz. 
In the baseline NASTRAN model, yrahjl70, the frequency is 22.08 Hz. This response is 
consistent with flight data and is characterized by vertical motion of the horizontal 
stabilizer and lateral response of the CPG station and pilot overhead. 
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Aft Fuselage Lateral Torsional Mode, 6.00 Hertz 




Figure 6.2: Aft Fuselage Lateral Torsion Mode Shape. 
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Horizontal Stab Yaw Mode, 1 1 .60 Hertz 
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Figure 6.3: Horizontal Stabilizer Yaw Mode Shape. 
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Horizontal Stab Vertical Bending Mode, 26.02 Hertz 






Figure 6.4: Horizontal Stabilizer Vertical Bending Mode. 
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C. FREQUENCY RESPONSE FUNCTIONS 



Frequency response functions were determined between the horizontal and vertical 
tail and accelerometers. To compute the FRFs, accelerometer locations had to be 
correlated to NASTRAN grids. Sikorsky provided a list, contained in Appendix D, of 
accelerometer locations and suggested NASTRAN grids. A partial list of the NASTRAN 
grid points used is shown in Table 6.2. 



Accelerometer/: 

Correlai 


NASTRAN 

tion 


Location 


NASTRAN Grid 


Nose Cone 


61641 


Pilot Heel 


6505 


Pilot Floor 


7246 


Copilot Heel 


7608 


Copilot Floor 


8948 


Pilot Overhead 


7514 


Copilot Overhead 


8614 


Vertical Pylon Top 


32341 


Hor Stab Left Tip 


41108 


Hor Stab Right Tip 


48108 



Table 6.2: Accelerometer Locations in NASTRAN Model. 

A one pound force was applied at these nodes with a frequency sweep of 5-25 Hz. 
Accelerations, in g’s and in all three directions, were found at each accelerometer location. 
The FRFs were written to a .pch file, converted into matrices compatible with MATLAB, 
and stored as .mat files for use in solving for PSDs and CSDs of loads. 



* 
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1 . 



Horizontal Tail 



The NASTRAN nodes along the quarter chord of the horizontal tail are shown in 
Figure 6.5. Only one loading condition, loads applied in the z direction, was used. In 
addition to determining the FRF at each grid point along the quarter chord and the 
accelerometer locations, a chordwise investigation was conducted. At five different span 
locations, the FRF was found for each node along the chord. This provided no additional 
information as the FRF 1 varied minimally from leading edge to trailing edge. Therefore, it 
is not included. 




Figure 6.5: Horizontal Tail Nodes and Loading Condition. 
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Figures 6.6 and 6.7 show the FRF between the horizontal stabilizer and the pilot 
heels lateral and vertical. Although the FRF was determined from 5-25 Hz, the magnitude 
beyond 20 Hz was extremely small. For better resolution in the plot, the axis was reduced 
to 5-20 Hz. The span axis represents every node along the quarter chord, with the right 
tip defined as 0% span and the left tip representing 100%. As expected, the peak FRFs 
are aligned with the natural frequencies. 

In Figure 6.6, the lateral FRF has the greatest magnitude at 6 Hz on the right tip of 
the stabilizer. All other FRFs are less than this value. As the load was moved along the 
tail, the magnitude continued to decrease until a minimum at mid-span. This is expected 
as loads directly over the vertical tail induce very small moments, and the lateral motion 
throughout the fuselage reduces. Continuing out the left horizontal, the magnitude 
continues to increase; however, never reaches the amplitude from the right side. 

Figure 6.7 shows the FRF between the applied load and the vertical acceleration at 
the pilot heels. The FRF does have symmetry across the tail; however, the magnitude is 
substantially less than the lateral. 

Figures 6.8 and 6.9 display the FRF between the same points but at the discrete 
frequencies of 6 and 1 1 .6 Hz. The magnitude was normalized to the maximum amplitude 
FRF. Comparisons can now be more easily drawn for relative values of each FRF. These 
plots indicate the significant role the right side of the horizontal tail plays in transferring 
energy at both 6 and 11.6 Hz. 
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Figure 6.6: FRF Horizontal Tail to Pilot Heel - Lateral Response. 
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x 1 0 ” Horizontal Tail - Pilot Heel, Vertical Response 
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Figure 6.7: FRF Horizontal Tail to Pilot Heel - Vertical Response. 
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Horizontal Tail (c/4) - Pilot Heels, Lateral 
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Figure 6.8: FRF Horizontal Tail to Pilot Heel - Lateral. 
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-o 





0% 12% 28% 46% 54% 64% 78% 94% 



% Span - Right Tip to Left Tip 



Figure 6.9: FRF Horizontal Tail to Pilot Heel - Vertical. 
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2 . 



Vertical Tail 



The loading condition and NAS TRAN nodes along the quarter chord of the 
vertical tail are shown in Figure 6. 10. Again, only one loading condition, loads applied in 
the y direction, was used. As with the horizontal tail, a chordwise investigation was 
conducted. Similarly, the FRF chordwise at five span locations varied only slightly. 

Those results are not presented. 



C/4\ 



Vertical Tail Loading Condition 
Lateral Force 
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* Figure 6. 10: Vertical Tail Nodes and Loading Condition. 
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Figures 6. 1 1 and 6. 12 show the FRF between the vertical stabilizer and the pilot 
heels lateral and vertical. As with the horizontal, no significant response level was evident 
in the 20-25 Hz range. The span axis represents every node along the quarter chord, with 
the base of the tail defined as 0% span and the top representing 100%. 

In Figure 6. 1 1, the lateral FRF of greatest magnitude is at 6 Hz. This occurs at the 
top of the vertical tail. All FRFs decreased as the load was moved down the tail, except 
at 11.6 Hz. The FRF for this frequency reached a maximum at the base of the vertical tail. 

Figure 6. 12 shows the FRF between the applied load and vertical acceleration at 
the pilot heels. It is consistent with the information from the previous plot. Magnitudes at 
all but 1 1 .6 Hz decrease as the load was moved down the tail. The amplitudes are 
substantially less than lateral. 

Figures 6.13 and 6.14 display the FRF between the same points but at the discrete 
frequencies of 6 and 1 1.6 Hz. As with the horizontal stabilizer, the magnitude was 
normalized to the FRF of maximum amplitude, which occurred on the right tip of the 
horizontal stabilizer. 

The top of the vertical has the same magnitude at 6 Hz as does the right tip of the 
horizontal. The base of the vertical tail has the dominant FRF at 11.6 Hz. Comparing this 
with the pressure transducer data, the right horizontal tail had the most energy 
concentrated in the 5-8 Hz bandwidth. This is where the 6 Hz FRF is a maximum. For 
the base of the vertical tail, the pressure energy was broadband from 11-19 Hz. The 11.6 
Hz FRF is a maximum at this location. 
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Figure 6. 1 1 : FRF Vertical Tail to Pilot Heel - Lateral Response. 
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g' s/lb 



3 Vertical Tail - Pilot Heel, Vertical Response 




Figure 6. 12: FRF Vertical Tail to Pilot Heel - Vertical Response. 
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Vertical Tail (c/4) - Pilot Heel, Lateral 




Figure 6. 13: FRF Vertical Tail to Pilot Heel - Lateral. 
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Figure 6. 14: FRF Vertical Tail to Pilot Heel -Vertical. 
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VII. RESULTS AND ANALYSIS 



A. SOLUTION METHOD 

Having obtained the PSDs of the response and the structural frequency response 
functions, a solution can be obtained for the spectral content of the applied loads. The 
method for accomplishing this is provided by Bendat and Piersol [Ref. 8:p. 240], A 
summary of this method follows. 

X is defined as a column vector representing the Fourier transform of q input 
records, and Y the Fourier transform of k output records, 

















y 2 


X = 


* 


Y = 


: 








Y k _ 



X*, Y* = complex conjugate of the vectors X, Y 
X', Y' = transpose of the vectors X, Y 

The input and output spectral density matrices can then be defined as 



G 

G 



xx 



yy 



jE^CX’} = input spectral density matrix 



Y'}= output spectral density matrix 
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The notation used by Bendat and Piersol for Gxx and Gyy simply means that the 
PSDs and CSDs of each input and output record are contained with these matrices. 
Another method for representing these terms, and useful for programming is 







G n 


G n •" 


G \ q 




G « 


- 


G 2, 


G 22 — 


G . 2 * 








G * 


G ,2 


G 99. 








G y,y, 


G y\yi 


... g 


y\yu 




- 


G w* 


G y^ 


... g 


yiy* 






_ G y t y, 


G y*y 2 


**• G y*y t . 



Both Gxx and Gyy are Hermitian matrices, with G y = Gj* for all i and j. Each 
element in both matrices represents either the PSD or CSD for the inputs and responses. 
When the subscripts are the same, as with Gn, the value is a real number and represents 
the PSD for the load at point one. The other elements of the matrix are complex, such as 
G12 which is the CSD between the loads at point one and two. The same is true for the 
response matrix, with the diagonal being the PSDs for responses and the off diagonal 
terms are the CSDs between responses. 

Finally, the FRF, Hxy, where x defines the input and y is the output is given by 
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Depending on the number of inputs and outputs used in the solution set, the FRF matrix 
may not be square. Using this definition of H X7 , it follows that the output Y is related to 
the input X by the relationship 

Y = H^,X where H^, is the transpose matrix to 

Having defined these terms, the following derivation is required to arrive at a 
solution for Gxx- 



y = h;x 

Y' = (H'^Xy = X'H^ 

Y* = (H^X)* = H^X* 

Y*Y = (H' xy *X*)(X'H xy ) 

Taking the expected value of both sides and multiplying by (2/T) gives 

g„=(h'„')g„(hJ 

It follows that 

Assuming the required inverse matrices exist, a solution can be obtained for the 
spectral content of the applied loads. In the case of unequal inputs and responses, where 
Hxy is not square, the solution, Gxx, is found using gaussian elimination. Conditioning of 
the matrices becomes critical during the computations. 
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B. MATLAB® PROGRAM 



Using this formulation, a MATLAB® code was written, contained in Appendix E, 
that solved for Gxx First, the user specified the frequency range of interest. This was 
limited from 5-25 Hz because of the range of frequency sweeps use to determine the 
frequency response functions. Then, NASTRAN grid points were selected for application 
of forces. Finally, response points to be used in the solution set were chosen. 

The program looks up, on a frequency by frequency basis, the stored FRFs 
between force application points and response locations. This created the Hxy matrix. 

The Gyy matrix was found, again on a frequency by frequency basis, from the selected 
response points. PSDs and CSDs for each set of accelerometer flight test data was stored 
in .mat files. This allowed retrieval for use in the appropriate Gyy matrix. Because the 
frequency interval on FRFs was .05 Hz and on the 1024 point FFT PSDs was .293 Hz, 
cubic spline interpolation was used to match the data sets. 

The solution, G^, represented the PSDs and CSDs of the applied loads at discrete 
frequencies. The frequency band was dependent on what the user had selected for the 
upper and lower limit. The interval between frequencies, or effective bandwidth was 
based on Gyy. Hence for this output, with a 1024 point FFT, the interval was .293 Hz. 
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C. RESULTS 



Numerous cases, with different locations of forces and reponses, were computed. 
Some generalities included all forces placed along the quarter chord of either the 
horizontal or vertical tail. These points are shown in Figures 6.5 and 6.10 and listed in 
Table 7. 1 . Additionally, the only response data used was from a 100 knot record, entitled 
run 24. The response points contained in this data set are also displayed in Table 7. 1 . 
Horizontal tail span is defined from the right tip. Vertical tail span begins at the base of 
the tail. 



Solution Parameters 


Response Data 


Force Application Points 


Run 24: Accelerometers 


Aircraft Location 


% Span 


NASTRAN 

Grid 


Nose Cone Lateral 


Horizontal Tail 
C/4 


12.38 


48008 


Nose Cone Vertical 


28.39 


46008 


Pilot Heel Lateral 


45.84 


45108 


Pilot Floor Lateral 


58.97 


44008 


Pilot Floor Vertical 


70.70 


43008 


Pilot Overhead Lateral 


86.70 


42008 


Copilot Overhead Lateral 


Vertical Tail 
C/4 


80.66 


31844 


Vertical Pylon Top Lateral 


65.05 


31842 


Hor Stab Left Tip Long 


49.45 


31840 


Hor Stab Right Tip Long 


33.85 


31838 


Hor Stab Left Tip Vertical 


18.25 


31836 


Hor Stab Right Tip Vertical 


6.81 


30722 



Table 7. 1 : Response Data and Force Application Points. 



The first solution set used two forces, one on the right horizontal stabilizer at mid 
semi-span and the other at mid-span of the vertical tail. Then, additional unknown forces 
were included, with different response data. For several frequency bands, the condition 
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number of the FRF matrix, Hxy, became exceedingly high. This occurred when the FRFs 
were very small or FRFs from different points were of equal amplitude. In these cases, the 
precision of MATLAB® could not distinguish between numbers and interpreted the two 
rows of the matrix as duplicates of each other. Therefore, the inverse was not reliable. 

Using between two and twelve forces, the response locations were varied to 
produced different solutions. Table 7.2 through 7.8 display the results for seven different 
loading conditions. In all these cases, every response location that run 24 contained was 
used. The tables present the RMS for each force within the frequency bands previously 
identified. In every case, the energy within these three bins represented over 80% of the 
total energy, throughout the 5-25 Hz range, of the force at that node. The final column 
indicates the RMS of the total force. 

The RMS values were determined for each frequency band and then divided by 
the largest RMS value of total force. They represent the percentage of the total force that 
is located at that point on the aircraft and within the frequency range. Similarly, the RMS 
of total force was normalized to this same value. Hence, the RMS of total force that 
equals 1.0 represents the location on the tail that had the largest force throughout the 5-25 
Hz band. Normalizing the forces was done to aid in comparison of loads between cases. 
Additionally, the conditioning of the matrices produced highly inaccurate results within 
certain frequency bands. In these cases, the solution had large amplitude PSDs in 
frequency ranges containing little response and small FRFs. The loads in nine and twelve 
force cases. Table 7.7 and 7.8, are highlighted to show the frequency content and 
magnitudes of the computed loads. 
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Forcing Function - Normalized (Ib/lb) 




Frequency Bin 


RMS of 
Total Force 


Node 


5-8 hz 


10-1 3 hz 


22-25 hz 


48008 


0.34 


0.38 


0.21 


0.44 


31840 


0.47 


1.00 


0.52 


1.00 



Table 7.2: RMS of Two Force Solution Set. 



Forcing Function - Normalized (Ib/lb) 




Frequency Bin 


RMS of 
Total Force 


Node 


5-8 hz 


10-1 3 hz 


22-25 hz 


42008 


0.36 


0.42 


0.11 


0.60 


46008 


0.39 


0.39 


0.20 


0.71 


31840 


0.46 


0.71 


0.30 


1.00 



Table 7.3: RMS of Three Force Solution Set. 



Forcing Function - Normalized (Ib/lb) 




Frequency Bin 


RMS of 


Node 


5-8 hz 


10-13 hz 


22-25 hz 


Total Force 


43008 


0.14 


0.06 


0.12 


0.22 


46008 


0.15 


0.19 


0.09 


0.28 


31844 


0.51 


0.49 


0.36 


0.88 


31836 


0.76 


0.38 


0.26 


1.00 


Table 7.4: RM 


S of Four Force Solution Set. 
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Forcing Function - Normalized (Ib/lb) 




Frequency Bins 


RMS of 
Total Force 


Node 


5-8 hz 


10-13 hz 


22-25 hz 


43008 


0.22 


0.15 


0.07 


0.30 


45108 


0.53 


0.28 


0.06 


0.63 


46008 


0.30 


0.16 


0.04 


0.36 


31844 


0.50 


0.16 


0.10 


0.60 


31836 


0.89 


0.21 


0.13 


1.00 



Table 7.5: RMS of Five Force Solution Set. 



Forcing Function - Normalized (Ib/lb) 




Frequency Bins 


RMS of 
Total Force 


Node 


5-8 hz 


10-13 hz 


22-25 hz 


42008 


0.03 


0.01 


0.00 


0.03 


44008 


0.07 


0.02 


0.01 


0.07 


46008 


0.03 


0.01 


0.00 


0.03 


31844 


0.48 


0.11 


0.05 


0.47 


31840 


1.00 


0.27 


0.09 


1.00 


31836 


0.53 


0.17 


0.05 


0.54 



Table 7.6: RMS of Six Force Solution Set. 



Forcing Function - Normalized (Ib/lb) 




Frequency Bin 


RMS of 
Total Force 


Node 


5-8 hz 


10-13 hz 


22-25 hz 


46008 


0.00199 


0.00564 


0.00034 


0.00613 


45108 


0.01034 


0.00795 


0.00135 


0.01426 


44008 


0.01581 


0.01946 


0.00173 


0.02596 


43008 


0.00661 


0.01301 


0.00036 


0.01471 


31844 


0.19509 


0.15047 


0.00907 


0.25202 


31842 


0.40794 


0.32484 


0.01723 


0.53204 


31838 


0.58917 


0.56503 


0.01648 


0.83418 


31836 


0.60257 


0.75681 


0.02316 


1.00000 


30722 


0.23072 


0.36292 


0.01422 


0.44914 



Table 7.7: RMS of Nine Force Solution Set. 
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Forcing Function - Normalized (Ib/lb) 




Frequency bin 


RMS Of 
Total Force 


Node 


5-8 hz 


10 - 1 3 hz 


22-25 hz 


48008 


0.00027 


0.00777 


0.00005 


0.00818 


46008 


0.00053 


0.02362 


0.00015 


0.02487 


45108 


0.00047 


0.02538 


0.00018 


0.02738 


44008 


0.00148 


0.01474 


0.00021 


0.01830 


43008 


0.00122 


0.01450 


0.00018 


0.01670 


42008 


0.00028 


0.00363 


0.00004 


0.00415 


31844 


0.03813 


0.09873 


0.00017 


0.10854 


31842 


0.20775 


0.47014 


0.00037 


0.53079 


31840 


0.43149 


0.86229 


0.00037 


1.00000 


31838 


0.40227 


0.72978 


0.00030 


0.86572 


31836 


0.17668 


0.29943 


0.00009 


0.36249 


30722 


0.03467 


0.07152 


o.oooTo 


0.08449 



Table 7.8: RMS of Twelve Force Solution Set. 



Figure 7.1 through 7.4 show the PSDs of loads for the twelve force solution set. 
The loads on the horizontal tail are within frequency bins which have little response. The 
magnitudes are much smaller than forces on the vertical and associated with poorly 
conditioned Hxy matrices. Under the PSD plots is the phase angle between point 31840, 
where the force of largest magnitude is located, and selected points. They were 
determined from CSDs which capture the phasing information. Bendat and Piersol [Ref. 
8:p. 125] provide a method of determining the phase angle from cross-spectral density 
functions. Although the plots are difficult to interpret, if viewed at just the modal 
frequencies earlier identified, 6.0 and 1 1.6 Hz, the phase relationship can be seen. This is 
available on the vertical tail where forces at these frequencies are located. The horizontal 
tail has no forces at these frequencies. 
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Normalized PSD - lb /hz 



xIO ' 4 



PSD of Forces - Right Horizontal Stab 




Phase Angie 31840 - 45108 




Frequency - Hertz 



Figure 7.1: PSD of Forces on the Right Horizontal Tail. 
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Normalized PSD - lb /hz 




5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 
Phase Angle 31840 - 43008 




Frequency - Hertz 



Figure 7.2: PSD of Forces on the Left Horizontal Tail. 



105 



degrees Normalized PSD - ttr/hz 



PSD of Forces - Vertical Stab (50%-100% Span) 
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Figure 7.3: PSD of Forces on the Vertical Tail - 50% to 100% Span. 
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Normalized PSD - lb /hz 



PSD of Forces - Vertical Stab (0%-50% Span) 
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D. ANALYSIS 



Beginning with the initial load set and continuing through twelve forces, the 
computations place the loads of largest magnitude on the vertical tail. Initially, the 10-13 
Hz frequency bin had the greatest energy. As more and more forces were added, the 
magnitude of forces in the 5-8 Hz range grew tremendously. This coincided with a poorly 
conditioned FRF (Hxy) matrix at 5 Hz. Even when eliminating the contribution to RMS 
from this frequency, forces from the base of the vertical tail to mid-span, at both 5-8 and 
10-13 Hz, play the dominant role in producing the response. 

In going from nine to twelve forces, the 10-13 Hz frequency bin, again, became the 
dominant source of energy. The problems with the conditioning of the FRF matrix, 
discussed above, were again present in the final solutions. The point where the largest 
forces were placed was grid 31840. Figures 7.5 and 7.6 display the distribution of forces 
along the horizontal and vertical tail. It can be seen that more energy was placed on the 
lower two-thirds of the vertical tail in both the 5-8 and 10-13 Hz frequency band. This is 
consistent with the results from the nine load case, shown in Table 7.7. 

The results also match the information from the analysis of the frequency response 
functions and pressure transducers. The solution places loads on the fuselage at points 
that have high magnitude FRFs at 11.6 Hz. Pressure transducers confirm broadband 
excitation, including this frequency, and of largest amplitude at the base of the vertical tail. 

Although the solution indicates a force on the vertical tail in the 5-8 Hz frequency 
bin, pressure transducers do not support this. The solution also has no forces in the 22-25 
Hz range. The lack of NAS TRAN modes in this bin has previously been discussed. 
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Horizontal Tail (c/4) - RMS of Forces 




12.39% 28.39% 45.84% 58.97% 70.69% 86.71% 

%Span -Right Tip to Left Tip 



Figure 7.5: Distribution of Forces, 12 Force Solution Set — Horizontal Tail. 
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Figure 7.6: Distribution of Forces, 12 Force Solution Set - Vertical Tail. 
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VIII. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

The goal of this work was to quantify the magnitude and location of aerodynamic 
loads contributing to the random excitation of the empennage. Then, follow-on efforts 
would have the proper focus and direction in modifying the structure. This was 
accomplished by first conducting an extensive analysis of flight test response data and 
pressure transducers. Due to the nature of the accelerometer time histories, spectral 
analysis was used. With this accomplished, the mode shapes and frequency response 
functions of the aircraft were determined. ANASTRAN model, modified to replicate the 
test flight aircraft, provided this information. Finally, a MATLAB® program was written 
to compute the spectral content of the aerodynamic loads producing the random 
excitation. 

Each portion of the analysis provided insight into the nature of the response, 
locations of excitations, and the method in which the structure transfers input to output. 
From all of these elements, combined, conclusions were drawn 

The dominant response is lateral and occurs in three frequency bands or bins. The 
first bin is 5-8 Hz, then 10-13 Hz, and finally 22-25 Hz. Each accelerometer has its 
maximum value within one of these bins. The top of the vertical pylon has maximum 
response at 6.00 Hz. For the nose cone, pilot’s station, and horizontal tail, 10-13 Hz 
excitation is the largest. Finally, the 22-25 Hz response is restricted primarily to the CPG 
station and vertical acceleration of the horizontal tail. 
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Within the first two bins, there are mode shapes that correspond with the motion 
seen in the output of the 22 accelerometers. The aft fuselage lateral torsion mode at 6.00 
Hz and horizontal stabilizer yaw mode at 1 1.6 were identified not only through 
NASTRAN but also during actual frequency response hover testing. No mode was 
identified or significant frequency response function found within the 22-25 Hz bin. The 
horizontal stabilizer vertical bending mode at 26.02 Hz matches the response seen in the 
22-25 Hz bin. The actual natural frequency of this mode was not identified during flight 
test; however, in previous NATRAN models it is at 22.06 Hz. 

From the forces computed through the MATLAB® program, good correlation 
existed with the previous analysis. Numerous loading combinations were applied with 
maximum loads continuously placed on the vertical stabilizer from the base to 50% span. 
The 10-13 Hz frequency bin contained the largest force. This agreed with data from 
differential pressure transducers which showed the highest broadband pressures at the base 
of the vertical tail. 

From this, it can be concluded that forces within the 10-13 Hz bandwidth 
impinging on the lower two-thirds of the vertical tail play the most significant role in 
creating the random excitation. This is supported by the differential pressure transducers. 
Solutions want to place 5-8 Hz forces in this same area. This is not supported by pressure 
transducers and leads to the conclusion that main rotor IP is the primary source of 6 Hz 
excitation. Because the 22-25 Hz excitation is so localized and no FRFs of magnitude are 
in this region, the least squares solution provides no contributing loads. If the horizontal 
stabilizer vertical bending mode is not at 26.02 Hz but lower, the broadband excitation 
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seen in the differential pressure transducers can excite this mode and produce the recorded 
response. 

B. RECOMMENDATIONS 

1. Correlate with Wind Tunnel Testing 

Boeing and Sikorsky have conducted two wind tunnel tests. The purpose was to 
quantify the flow impinging on the tail and identify changes to the surface of the aircraft, 
such as strakes, that would modify the flow and reduce excitation levels. The results of 
this analysis should be compared to data from the wind tunnel. This could confirm the 
conclusions of this study about the nature of the aerodynamic loads and their primary 
locations on the fuselage. 

2. Parameter Identification and Correlation with NASTRAN 

During hover frequency response testing, actual natural frequencies were identified 
for three modes. As additional test flight hours become available, more parameter 
identification can occur. This will better identify the exact frequency of the horizontal 
stabilizer vertical bending mode, as well as clarify other modes. In addition, the change in 
natural frequencies due to adding stiffening brackets and tip weights to the horizontal 
stabilizer can be assessed. With this information, a detailed analysis of the NASTRAN 
model can be conducted with modification of natural frequencies and mode shapes to 
replicate the actual aircraft. 
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3. Additional Airspeed Analysis 

The computed solution using the method of Bendat and Piersol was done for only 
100 knot records. After receiving additional airspeed records from Sikorsky, the PSDs 
and CSDs were computed; however, no complete solutions were obtained for the spectral 
content of forces. Calculating solutions throughout the range of airspeeds can provide 
further confirm and clarity as to when the random excitation begins and if the 6 Hz 
response is due to main rotor IP. 

4. Structural Optimization 

Ultimately, the goal is to eliminate or reduce the response to acceptable levels. 
Using the conclusions from this analysis and confirmation from the sources discussed 
above, the structure can be altered to reduce the frequency response functions within the 
excitation bandwidth. This can be accomplished using several different methods, or, 
optimally, a combination of all. 

First, using the NASTRAN model, structural modifications can be done that move 
the natural frequency, 11.6 Hz, of the horizontal stabilizer yaw mode. The focus of the 
changes should occur on the lower two-thirds of the vertical tail. Extensive work was 
done by Tobin [Ref. 12] and Shoop [Ref. 13] to increase torsional stiffness of the aft 
tailcone. The work attempted to stiffen the tail and move the aft fuselage lateral torsional 
mode away from MR IP. Their analysis and strain energy diagrams provide insight to 
methods of altering the structure, specifically for the 1 1 .6 Hz mode. 

Second, changes are proposed to the design of the tail fold, spars of the vertical 
tail, and the tail fan gear box struts. Additionally, a trade study is being conducted by the 
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Comanche Program Management Office (PMO - Comanche) on the mounting of the tail 
landing gear bay. As these designs are finalized, the dynamics of the tail will change. The 
quantification of the flow determined through this analysis should aid in creating a 
dynamically “tuned” design. 

Finally, modifying the flow around the main rotor hub and pylon has been 
investigated through the wind tunnel experiments previously discussed. Modifications to 
the structure must be done in conjunction with changes to the flow. Only through this 
integrated approach can an acceptable solution be obtained. In this manner, both the flow 
and structural response can be optimized for maximum reduction of the airframe vibratory 
response. 
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APPENDIX A: TEST FLIGHT PARAMETERS 



... 

Mnemonic 


Parameter 


Affected Components 


Units 


HBNORML1 


Horz Stabilizer Flatwise 
Bending, LHS (near Root) 


Horizontal Stabilizer 


in-lb 


HBNORMR1 


Horz Stabilizer Flatwise 
Bending, RHS (near Root) 


Horizontal Stabilizer 


in-lb 


HBNORML2 


Horz Stabilizer Flatwise 
Bending, LHS (outboard) 


Horizontal Stabilizer 


N/A 


HBNORMR2 


Horz Stabilizer Flatwise 
Bending, RHS (outboard) 


Horizontal Stabilizer 


N/A 


HSTBFWDB 


Horz Stab Fwd Folding Joint 
Bending LHS 


Horizontal Stabilizer 


in-lb 


HSVAFTL 


Horz Stab Aft LHS Spar Shear 
Force - lb 


Horizontal Stabilizer 


N/A 


HSVAFTR 


Horz Stab Aft RHS Spar Shear 
Force - lb 


Horizontal Stabilizer 


N/A 


HSVFWDL 


Horz Stab Forward LHS Spar 
Shear Force - lb 


Horizontal Stabilizer 


N/A 


HSVFWDR 


Horz Stab Forward RHS Spar 
Shear Force - lb 


Horizontal Stabilizer 


N/A 


VSFFTS1 


Vert Stab Attach Fwd Ftg Outer 
Cap LHS 


Vert / Shroud 


u‘7” 


VSFFTS2 


Vert Stab Attach Fwd Ftg Outer 
Cap RHS 


Vert / Shroud 


u‘T 


VSAFTS1 


Vert Stab Attach Aft Ftg Outer 
Cap LHS 


Vert / Shroud 


u‘7” 


VSAFTS2 


Vert Stab Attach Aft Ftg Outer 
Cap RHS 


Vert / Shroud 


u l T 



Table A. 1 : Test Flight Parameters. 
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V- 

Mnemonic 


.w.:- • ■ : ' 

Parameter 


Affected Components 




BVSRTMX 


Vertical Stab Root Roll 
Moment 


Vert / Shroud 


N/A 


DHSDRAG 


Horizontal Stab Drag - DP 


Horizontal Stab 


N/A 


DHSTABMX 


Net Stab Roll Moment -DP 


Horizontal Stab 


in-lb 


DHYAWMZ 


Horizontal to Vertical Stab 
Yaw Moment - DP 


Horizontal Stab 


in-lb 


DHSTABVZ 


Net Stabilizer Lift (lb) - DP 


Horizontal Stab 


N/A 


DVSAFTLH 


Derived Vert Load Aft LHS 
Lug - DP 


Vert / Shroud 


lb 


DVSAFTRH 


Derived Vert Load Aft RHS 
Lug - DP 


Vert / Shroud 


lb 


DVSFWDLH 


Derived Vert Load FWD LHS 
Lug - DP 


Vert / Shroud 


lb 


DVSFWDRH 


Derived Vert Load FWD LHS 
Lug - DP 


Vert / Shroud 


lb 


DVSRTMX 


Vertical Stab Root Roll MMT 
-DP 


Vertical Stab 
Attachment 


N/A 


DVSRTMY 


Vertical Stab Root Pitch 
MMT -DP 


Vertical Stab 
Attachment 


N/A 


DVSTABVY 


Derived Vertical Stab Lateral 
Force - DP 


Vert / Shroud 


N/A 


HBEDGEL1 


Horz Stab Edgewise Bending, 
LHS (near root) 


Horizontal Stab 


N/A 


HBEDGER1 


Horz Stab Edgewise Bending, 
RHS (near root) 


Horizontal Stab 


N/A 



Table A.2: Test Flight Parameters. 



* 
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APPENDIX B: PSD AND CSD PROGRAMS 



% psd_flt_data .m 

% Program reads the input flight data for a run and the computes the 
% PSD for each accelerometer time history 

% The input record is an 18,000 X n matrix with each column 
% representing another transducer 

% Input the data set you want to read, sampling rate, data set you 
% want plotted 

clear 

run_labels__file 

dt = 1/300; 
string =[ , runl40 1 ]; 
speed =[’140 Knots’]; 
eval ( [ ' load ' , string, 1 ; ' ] ) ; 
i = 0; 

for col = 4:32 
i = i + 1; 
n = 1024; 

noverlap = round (. 67*n) ; 

eval(['x = ', string, '_data {:, col) ;*]) ; 

% Using the SPEC command 

% [ Pxx, F] = spec (x, n, noverlap, hanning (n) , 1/dt) ; 

% Using the PSD command 

[Pxx, F] = psd (x, n, 1/dt , hanning (n) , noverlap) ; 

% Normalize by multiplying by 2/frequency sampling rate 
% To get the rms multiply by sampling rate/nfft then add 
% This is from Pareseval's theorem 

Pxx = (2/ (1/dt)) * Pxx; 

% Another way to normalize so that the rms is just amplitudes add 
% together 

% With this method, the value from different length fft's is the same 
% Pxx = (2/ (1/dt)) * ( 300/n) * Pxx; 
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% Plotting Function 



plot_max = round (25/ (( 1/dt) /n) ) ; 

ymax = max ( Pxx ( 10 : plot_max, 1 ) ) ; 

figure (i) 

plot ( F, Pxx ( : , 1 ) ) ; 

grid; 

axis([0 25 0 ymax]); 

set (gca, ' xtick 1 , [0:1:25]) ; 

eval ( [ ' title ( 1 1 PSD 1 , deblank (run_labels (col, : ) ) , 1 ' , num2str (n) , ' Point 

FFT 1 ! ) ' ] ) ; 

xlabel ( 1 Frequency - Hz'); 
ylabel ( ! g's^2/Hz 1 ) ; 

eval ( [ 1 text ( 5, . 9*ymax, * 1 Run 1 , speed, . . . 

1 , Hanning Window, 67% Overlap ! 1 ) 1 ]) ; 

% Save the PSD data in another matrix for use in solving for forcing 
functions 

% Name of the matrix saved is psd_run24 (mat format) and has the 
frequency and 

% the psd for all accelerometer data we had for run24. 



eval ( [ 1 psd_ f , string, 1 ( : , i) = Pxx (:,!);*]); 
end 

%save frequency as a separate variable 
psd_labels = run_labels ( 4 : size ( run_labels, 1 ) , : ) ; 

eval([ T save psd_ 1 , string, 1 psd_ string, 1 F psd_labels ' ] ) ; 



% 
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% csd_f 1 t_data . m 

% Program reads the input flight data for a run and the computes the 
% CSD (cross spectral density) between output accelerometer time 
% history 

% The input record is an 18,000 X n matrix with each column 
% representing another transducer 

% Input the data set you want to read, sampling rate, data set you 
% want plotted 

dt = 1/300; 
load run24 
k = 0; 

for col = 5:5 
for i = 1:12 
k = k + 1; 
ii = i + 3; 
n = 1024; 

noverlap = round (. 67*n) ; 
x = run24_data ( : , col) ; 
y = run24_data ( : , ii ) ; 



% Using the CSD command 

[Pxy, F] = csd (x, y,n, 1/dt, hanning (n) , noverlap) ; 

% Normalize by multiplying by 2/frequency sampling rate 
% To get the rms multiply by sampling rate/nfft then add 
% This is from Pareseval's theorem 

Pxy = (2/ (1/dt)) * Pxy; 

% Another way to normalize so that the rms is just amplitudes added 
% together 

% With this method, the value from different length fft's is the same 

%Pxy = (2/ (1/dt))* (300/ n) * Pxy; 

% Plotting Function 

plot_min = round (5/ (( 1/dt ) /n) ) ; 

plotjmax = round(25/ ( (1/dt) /n) ) ; 

ymax = max (abs ( Pxy (plot_min: plotjmax, 1 ) ) ) ; 

figure (k) 

subplot (2,1,1) ; 

plot (F, abs (Pxy ( : , 1) ) ) ; 

grid; 

axis ([5 25 0 ymax]); 

set (gca, ' xtick * , [5:1:25] ) ; 

eval ([ 1 title ('* CSD ', deblank (run24_labels (col, :)), ' - ', 
deblank ( run24_labels (ii, : ) ) , 1 1 1 ) 1 ] ) ; 

%xlabel (' Frequency - Hz'); 
ylabel ( 1 g' s A 2/Hz (Amplitude) 1 ) ; 
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subplot (2,1,2) 

plot ( F, 180/pi * angle ( Pxy (:, 1) )) ; 

angle_max = max (18 0/pi * angle ( Pxy (plot_min: plot_max, 1 ))) ; 

axis ( [5 25 -180 180] ) ; 

set (gca, 1 xtick ' , [5:1:25] ) ; 

set (gca, 1 ytick * , [-180:60:180] ) ; 

grid; 

xlabel (' Frequency - Hz * ) ; 
ylabel (* Phase Angle - degrees 1 ) 

eval ( [ 1 title ( 1 1 Run ' , num2str ( run24_data (1,1)),... 

* , Hanning Window, 67% Overlap, * , num2str (n) , * Point FFT 1 ' ) 1 ] ) ; 
%eval ( [ ' text (5, . 89*angle_max, 1 1 Phase Angle 
* ,num2str (run24_data (1, 1) ) , . . . 

% * , Hanning Window, 67% Overlap, ' , num2str (n) , 1 Point FFT 1 * ) 1 ] ) 

% Save the CSD data in another matrix for use in solving for forcing 
% functions 

% Name of the matrix saved is psd_run24 (mat format) and has the 
% frequency and the psd for all accelerometer data we had for run24. 

eval ( [ 1 csd_run24_* ,num2str (col-3) , * ( : ,i) = Pxy( :,!);*]); 

end 



%eval ( [ 'save psd_run24 csd_run24_* , num2str ( col-3 ), ' F -append;*]); 
end 



% 
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APPENDIX C: MODE SHAPES 



Mode shapes for selected modes are presented on the next eleven pages. 
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Fuselage First Vertical Bending Mode, 7.70 Hertz 








; 1111 












HI -iilli 


11111 illlt ■ 


|| 





^ in 
o> co 
°o co 

co co 



CD 

h- 

CO 

CO 



in 



co 

o> 



CO 



co in 
c\i c\i 



o> 

CO 
CO 
c\j c\i 



o 

oo 

o 



CM 

00 



cm 

CD 

in 



CO 

o 

CO 



o 



00 1^ 
Tt in 

00 CNJ 

r- in 



CD 
CD 

cd in 

CM h- 
O 
O 








■ 


§|j| 


ill 




ilHij 


i|| 




|| 


;:;-v --l 


1 


I 


[ 


J|| 






in 


CD 




I s - 


00 


o> 


o 




CM 


CO 




00 


N- 


CD 




o> 


co 


h- 




in 


o> 


CO 


00 


CM 


CD 


o 




''tf' 


in 


CD 




00 


CD 


CO 




co 


in 


CO 


o 


00 


in 


CO 


o 


00 


CM 


CD 


in 


CO 


CO 


CO 


CO 


CM 


CM 


CM 


CM 


i- 




1- 




h- 


in 


CM 






nI ^ 






Horizontal Stab Roll Mode, 9.50 Hertz 
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Horizontal Stab Yaw Mode, 1 1 .60 Hertz 
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2nd Fuselage Vertical Bending Mode, 14.52 Hertz 
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Fuselage Mode, 20.45 Hertz 
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Stab Vertical Bending Mode, 21 .56 Hertz 
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APPENDIX D: ACCELEROMETER LOCATIONS & NASTRAN GRIDS 



RAH66- MODEL J yrahj194.dat 4/17/98 sp 

COMPARISON OF KEY FLIGHT TEST LOCATIONS WITH SUGGESTED NASTRAN APPROXIMATE LOCATIONS 



ITEM 


LOCATION 


Mnemonic 


FLIGHT TEST 
Approximate Location 
STA BL WL 

mm mm mm 


NASTRAN 
Location 
STA BL 

mm mm 


WL 

mm 


NASTRAN 

GRID 

POINT 


1 


NOSE 


VNOSE:D:300 


5286.00 


0 00 


2720.00 


5298.00 


0.00 


2720.00 


61641 






LNOSE:D:300 




















LONOSE:D:300 
















2 


PILOT HEEL 


LPHEELD300 


6452.00 


150.00 


2720.00 


6306.00 


300.00 


2720.00 


6205 






VPHEELD300 








6306.00 


122.00 


2720.00 


6245 














6692.00 


300.00 


2720.00 


6505 














6692.00 


122.00 


2720.00 


6545 


3 


PILOT FLOOR 


LPF:D300 


7450.00 


000 


2720.00 


7525.00 


0.00 


2720.00 


7600 






LOPF.D300 








7340.00 


122.00 


2720.00 


7245 






VPF:D300 








7340.00 


-122.00 


2720.00 


7246 














7525.00 


110.00 


2720.00 


7645 














7525.00 


-110.00 


2720.00 


7646 


4 


COPILOT HEEL 


VCOPHEEL:D:300 


7740.00 


150.00 


2980.00 


7557.00 


110.00 


2838.00 


7659 














7589.00 


300.00 


2956.00 


7607 














7589.00 


-300.00 


2956.00 


7608 














8084.00 


300.00 


3025.00 


7907 














8084.00 


-300.00 


3025.00 


7908 


5 


COPILOT FLOOR 


VCOPF:D:300 


8900.00 


0.00 


3130.00 


9050.00 


0.00 


3160.00 


9100 














8871.00 


140.00 


3135.00 


8947 














8871.00 


-140.00 


3135.00 


8948 














8960.00 


140.00 


3147.00 


9047 














8960.00 


-140.00 


3147.00 


9048 


6 


PILOT OVERHEAD 


LPOVRHD:D:300 


7420.00 


0.00 


4100.00 


7349.00 


176.00 


4153.00 


7513 














7349.00 


-176.00 


4153.00 


7514 














7733.00 


227.00 


4229.00 


7713 














7733.00 


-227.00 


4229.00 


7714 


7 


COPILOT OVERHEAD 


LCPOVRHD:D:300 


8795.00 


200 00 


4200.00 


8602.00 


250.00 


4323.00 


8613 














8602.00 


-250.00 


4323.00 


8614 














8849 00 


251.00 


434300 


8713 














8849.00 


-251.00 


4343.00 


8714 


8 


VERTICAL PYLON TOP 


LVPYLNTD300 


17737.00 


185.00 


5360.00 


17737.10 


151.15 


5360.00 


32338 














17737.10 


159.90 


5360.00 


32341 














17743.30 


162.40 


5367.97 


91225 














17715.30 


164.67 


5370.54 


91236 


9 


HOR. STABILIZER LEFT TIP 


VSTABTLT:D:300 


17670.00 


-1100.00 


5463.00 


17670.10 


-1219.20 


5463.39 


41008 






LOSTABTL:D:300 








17774.40 


-1219.20 5463.20 


41007 














17670.10 -1068.90 5463.39 


41108 














17774.40 


-1068.90 5463.20 


41107 


10 


HOR. STABILIZER RIGHT TIP 


VSTABTRT:D:300 


17670.00 


1580.00 


5463.00 


17670.10 


1400.77 


5463.39 


49008 






LOSTABTR:D:300 








17774.40 


1400.77 


5463.20 


49007 














17666.70 


1400.77 


5374.97 


49001 














17770.70 


1400.77 


5364.18 


49002 


11 


HOR. STABILIZER VERTICAL 


HBNORML1 :D:300 




-275 00 




17670.10 


-281.41 


5463.39 


43108 




BEND. LT SIDE 










17666.70 


-281.41 


5374.97 


43101 














17982.30 


-281.41 


5444.27 


43105 














17979.90 


-281 .41 


5367.95 


43104 



Table D. 1 : Accelerometer - NASTRAN Locations. 
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ITEM 

No. 


LOCATION 


Mnemonic 


FLIGHT TEST 
Approximate Location 
STA BL WL 

mm mm mm 


NASTRAN 

Location 

STA BL WL 

mm mm mm 


NASTRAN 

GRID 

POINT 


12 


HOR. STABILIZER VERTICAL 


HBNORML2:D:300 


-660.00 


17670.10 -648.36 5463.39 


42108 




BEND. LT SIDE 






17666 70 -648.36 5374.97 


42101 










17982.30 -648.36 5444.27 


42105 










17979.90 -648.36 5367.95 


42104 


13 


HOR. STABILIZER VERTICAL 


HBNORMR1 :D:300 


380.00 


17670.10 341.60 5463.39 


45308 




BEND. RT SIDE 






17666.70 341.60 5375.00 


45301 










17982.40 415.31 5444.25 


45305 










17980.00 415.31 5367.95 


45304 


14 


HOR. STABILIZER VERTICAL 


HBNORMR2:D:300 


770.00 


17670.10 656.83 5463.39 


46008 




BEND. LT SIDE 






17666.70 656.83 5374.97 


46001 










17982.60 656.83 5444.24 


46005 










17979.90 656.83 5367.95 


46004 


15 


HOR. STABILIZER EDGE 


HBEDGEL1 :D:300 


-240.00 


17670.10 -281.41 5463.39 


43108 




BEND. LT SIDE 






17666.70 -281.41 5374.97 


43101 










17982.30 -281.41 5444.27 


43105 










17979.90 -281.41 5367.95 


43104 


16 


HOR. STABILIZER EDGE 


HBEDGER1 ;D:300 


350.00 


17670.10 341.60 5463.39 


45308 




BEND. RT SIDE 






17666.70 341.60 5375.00 


45301 










1798240 415.31 5444.25 


45305 










17980.00 415.31 5367.95 


45304 


17 


LT AFT SHROUD/V.TAIL 


VSAFTS1:D:300 




17002.00 -294.50 4400.00 


21326 




INTERFACE FTT1NG 










18 


RT AFT SHROUD/V.TAIL 


VSAFTS2:D:300 




17002.00 -176.00 4400.00 


21327 




INTERFACE FTTING 










19 


LT FWD SHROUD/V.TAIL 


VSFFTS1:D:300 




16635.00 -294.50 4400.00 


24621 




INTERFACE FTTING 










20 


LT AFT SHROUD/V.TAIL 


VSFFTS2:D:300 




16635.00 -176.00 4400.00 


24620 




INTERFACE FTTING 










21 


FWD H. STAB/V.TAIL • 


DHSTABMXD300 




17730.00 62.85 5404.00 


32391 




FrrriNG cl 


DHSTABMZ:D300 




17730.00 151.15 5404.00 


32394 






DHSTABVZ:D300 




17737.10 104.55 5360.00 


32336 



% 



Table D.2: Accelerometer - NASTRAN Locations. 
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APPENDIX E: LOAD CALCULATIONS PROGRAM 

% f f_solve_nonsym_inter .m 

% Program that solves for the forcing functions, in PSD by using the 

% techniques discussed by Bendat and Pierson for multiple input 

% single output 

% Gxx = (Hxy* 1 ) *-1 * Gyy * (Hxy)^-l 

% First must identify the location (nodes) of the forcing function 

% the location of the response (nodes), and the direction of the 

% response 

% The directions are somewhat backwards - due to the way I had Nastran 
% write out the punch files ( which became .mat files) 

% x = 3 (longitudinal) 

% y = 2 (lateral) 

% z = 1 (vertical) 

% The response data uses the following script: 

% First is node in Nastran Mode.l, next is direction and last is the 
% pointer into the psd matrix (which column) 



ncl = 


[61641 


2 


l] ; 


% 


Nose Cone Lateral 


phi = 


[6505 


2 


2] ; 


% 


Pilot Heel Lateral 


pfl = 


[7246 


2 


3] ; 


% 


Pilot Floor Lateral 


pol = 


[7514 


2 


4] ; 


% 


Pilot Overhead Lat 


cpol = 


[8614 


2 


5] ; 


% 


CPG Overhead Lat 


ncv = 


[61641 


1 


6] ; 


% 


Nose Cone Vertical 


pfv = 


[7246 


1 


7] ; 


% 


Pit Floor Vertical 


vptl = 


[32341 


2 


8] ; 


% 


Vert Pylon Top Lat 


hrtv = 


[48108 


1 


9] ; 


% 


HorStab Rt Tip Vert 


hltv = 


[41108 


1 


10] ; 


% 


HorStab Lt Tip Vert 


hrtlo = 


[48108 


3 


li] ; 


% 


HorStab Rt Tip Long 


hltlo = 


[41108 


3 


12] ; 


% 


HorStab Lt Tip Long 



% Use the following format to enter the input-ouput locations and 
% specify which run set 

run = 24; 

%no_fft = 4096; 

% force = [32311 46008 ] 

% response = [ pol 
% phi]; 

force = [ 42008 44008 46008 31844 31840 31836 ]; 
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response = [ phi 
pfv 
cpol 
ncl 
ncv 
vptl 
hrtv 
hltv 
hrtlo 
hltlo ] ; 

% Assumptions is that loads are correlation - are the loads correlated? 

% If the loads are completely uncorrelated then set corr = 0 
% If it is something inbetween then define the matrix relate 

force_corr = 1; 
response_corr = 1; 

% Use the following format for relate - forcing function (load) to load 
% (which are correlated) 

% force_relate = [110 
% 110 

% 001] This tells us force 1 is related to force 2, 

% all others are not 

% force_relate = [10 

% 0 1 ]; 

% Load the frequency response function matrix 
% Counter to use for matrix index 

k = 0; 



[no_row / force__num] = size (force); 

[ r espons e_num, no_col ] = size ( response) ; 



for f = l:force_num 

% frf file must have the name format frfnode 

for r = 1 : response_num 
k = k+1; 

eval ( [ 1 load f rf 1 , num2str ( force ( 1, f ) ) , ' ; * ] ) 

eval ( [ 1 frf_sel = f rf_ f , num2str ( response ( r, 1 )), T ; *]) 

y = [ (frf_sel(:,l) + frf_sel(:,2) .*i) ( frf_sel ( : , 3) + 

frf_sel(:,4) . *i) ( f rf_sel ( : , 5) + frf_sel(:,6) . *i)]; 

H(: f k) = y (:, response (r, 2) ) ; 
end ^ 



end 
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% H (frf matrix has now been loaded in the following format; from our 
% example nodes - 

% H = [ Hll H12 H13 H21 H22 H23 H31 H32 H33] ; 

% 

% Because of the weird way I saved the pch files, the vector {time} is 
% actually frequency 

% Now we need to load the response data from the PSD and CSD we have 
% done before 



eval (['load psd_run 1 , num2str ( run) , 1 ; ' ]) 

if exist ('no_fft') == 1 

eval ( [ ' load psd_run' , num2str ( run) , , num2str (no_fft) , 1 ; ' ] ) 

end 



k = 0; 

for f = 1 : response_num 

for r = 1 : response_num 
k = k+1; 
if f — r 

G(:,k) = psd_run24 (:, response (r, 3) ) ; 
else 

eval ( [ 1 G ( : , k) = ... 

csd_run24_* , num2str ( response ( f , 3) ) , * ( : , response (r, 3) ) ; 1 ] ) 

end 

end 

end 

% Form the Hxy and Gyy matrix for each individual frequency, 

% solve for Gxx 
% and save in another matrix 

kk = 0; 

start_freq = 5.50; 
end_freq = 24.00; 

[ f_dif f , f_index_start] = min ( abs ( s tart_f req - F) ) ; 

[ f dif f , f_index_end] = min (abs (end_f req - F) ) ; 



for ff = 


f index start : 


f index_ 


G count = 
k = 0; 


0; 




kk = kk + 


l; 




for 


1 = 1: response 


num 


for 


m = 1 : response 


num 



139 



oP 6IP 6P dP 



% Load the Gyy square matrix 

Gyy(l,m) = G (f f , G_count+m) ; 
end 

G_count = G_count + response_num; 
end 

for 1 = l:force_num 
for m = 1 : response^num 

Load the Hxy square matrix 

[min_dif , val_index] = min (abs ( F ( f f ) - time’)); 
Cubic Spline Interpolation 

Hxy (l,m) = interpl (time, H ( : , k+m) , F ( f f ) , 9 spline ’ ) ; 
% Linear Interpolation 

Hxy (1 ,m) = interpl (time, H ( : , k+m) , F ( f f ) , 1 spline ' ) ; 
end 



k = k + response_num; 
end 

if exist (* force_relate 1 ) == 1 
Hxy = force_relate . * Hxy; 
end 

if exist ( 1 response_relate 1 ) == 1 
Gyy = excite_relate . * Gyy; 
end 

if response_corr == 0 

Gyy = eye (size (Gyy) ) . * Gyy; 

end 



% Now we have to solve for the Gxx matrix via the technique 
% discussed above 

if force_corr == 0 & response__corr == 0 
Gyy = diag(Gyy); 

Gxx = ( (abs (Hxy 1 ) ) ."2) \ Gyy; 
else 

Gxx = Hxy'\ Gyy /Hxy; 
end 

% Saving the output data Gxx to a matrix of the following form 
% Gxx = [Gil G12 G13 G21 G22 G23 G31 G32 G33] 

% Also outputs a vector, Freq, that has the frequencies 
% each row represents 

[Gxx_row, Gxx__col] = size (Gxx); 

% 

if Gxx_col == 1 
G_output (kk, : ) = Gxx 1 ; 

else 

k = 0; 



140 



for 11 = l:Gxx_row 
for mm = 1 : Gxx_col 

G_output (kk, k+mm) = Gxx(ll,mm); 
end 

k = k+ force_num; 

end 

end 

frequency ( kk) = F(ff); 
cond_Hxy (kk) = cond(Hxy); 
rank_Hxy(kk) = rank(Hxy); 
end 

big_Kahuna = [frequency 1 G_output cond_Hxy' rank_Hxy']; 
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